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1  INTRODUCTION 


Currently,  x-ray  mammography  is  the  only  imaging  modality  used  in  clinics  for  detection  of 
microcalcifications.  This  important  tool,  however,  is  ionizing  and  its  sensitivity  is  greatly  reduced  in 
dense  breasts.  The  vision  of  this  research  is  to  develop  a  non  invasive  ultrasound-acoustic  method 
that  can  detect  microcalcification  in  breast  with  sufficient  accuracy.  The  general  goal  of  this  research 
is  to  build  an  experimental  system  based  on  Ultrasound  Stimulated  Acoustic  Emission  (USAE)  for 
in  vitro  imaging  of  breast  tissue  specimens  and  test  its  performance  in  detecting  microcalcification. 
The  first  task  of  this  research,  which  covers  the  first  year  and  a  part  of  the  second  year,  focuses  on 
system  development  and  optimization  for  detection  of  beast  microcalcification.  This  task  includes 
the  development  of  both  the  theoretical  model  and  experimental  setup  required  for  the  method.  This 
report  addresses  the  activities  performed  in  the  first  year  in  both  areas.  In  particular,  the  report 
includes  the  development  of  the  theoretical  model  representing  the  relationship  between  the  incident 
ultrasonic  wave  and  the  resulting  acoustic  emission  field,  improvement  of  the  laboratory  setup  to 
reduce  the  acoustic  noise  interference,  and  improvement  of  the  scanning  system. 


2  BODY 

In  this  section  we  describe  research  accomplishments  associated  with  task  1.  This  task  focuses 
on  two  general  areas:  (1)  Theory  and  model  development;  (2)  System  improvement.  Both  areas  have 
been  addressed  here.  In  our  work  on  system  improvement  we  often  refer  to  some  physical  phenomena 
and  their  effect  of  image  quality.  These  phenomena  are  briefly  described  in  the  following  table.  Later 
we  will  present  the  methods  we  have  used  to  counter  the  unwanted  effects  of  these  phenomena  on 
images. 

Table  l  Major  problems  in  the  experimental  setup  of  USAE  system 


Cause 

Ultrasound  standing  wave  effect 
Acoustic  noise 

Mutiple  reflection  of  acoustic  emission 
Multipath 


Effects  on  image 

Contrast  dependency  on  object-transducer 
distance 

Limited  image  contrast,  low  detectability 
Contrast  dependency  on  water  level  in  the  tank 
False  shades 


2.a  Transducer  Design 

2.a.l)  High-frequency  Transducer 

We  have  designed  a  9  MHz  confocal  transducer.  This  transducer  should  provide  better  spatial 
resolution  (about  200  microns)  than  our  3  MHz  transducer.  This  transducer  is  under  construction. 
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2.a.2)  Composite  Transducer 

In  order  to  detect  small  objects  in  a  viscoelastic  background,  we  need  to  reduce  the  background 
noise  in  the  medium  (water  tank  in  our  experiments).  We  have  found  that  one  such  noise  source  is  the 
acoustic  field  resulting  from  the  vibration  of  the  transducer  case  at  the  acoustic  emission  frequency 
A u).  The  confocal  transducer  used  in  our  experiments  is  comprised  of  a  center  disc  element  and  an 
outer  ring  both  constructed  on  a  single  piezoelectric  ceramic  cap.  A  possible  source  of  the  vibration 
is  coupling  between  the  elements.  An  alternative  piezoelectric  material  that  could  be  used  for  a 
transducer  with  lower  inter-element  coupling  is  a  piezoelectric  composite  material.  Based  on  this 
strategy,  we  designed  and  constructed  a  confocal  composite  transducer,  with  the  center  frequency, 
diameter,  and  focal  length  of  3.5  MHz,  25  mm,  and  7  cm,  respectively.  The  new  transducer,  although 
less  noisy,  unfortunately  cannot  be  used  at  the  same  power  level  as  the  ceramic  transducer  because 
of  its  limited  ability  to  dissipate  heat. 

2.a.3)  Annular  Array  Design 

We  have  noted  that  when  we  use  a  confocal  transducer  to  focus  two  CW  ultrasound  beams  on 
an  object  with  a  high  reflection  coefficient,  the  ultrasound  beam  reflects  back  and  forth  between  the 
transducer  and  the  object,  producing  a  standing  wave  pattern.  Multiple  reflections  can  be  strong 
because  the  object  is  located  at  the  focal  point  of  the  concave  transducer,  allowing  not  only  the 
incident  beam,  but  all  the  multiple  reflected  beams  to  be  focused  on  the  object.  Because  of  the 
standing  wave  effect,  the  acoustic  emission  field  of  the  object  becomes  excessively  dependent  on  the 
distance  between  the  object  and  transducer. 

To  reduce  the  standing  wave  effect  we  have  designed  a  multi-element  planar  annular  array 
transducer.  The  elements  are  divided  into  two  groups  to  produce  a  center  beam  and  an  outer  beam 
similar  to  those  of  our  conventional  (concave)  two-element  confocal  transducer.  Here  the  focusing 
is  achieved  by  introducing  the  proper  phase  to  the  elements  of  the  array  in  each  group.  Hence  in 
transmit  mode,  the  transducer  operates  similar  to  our  conventional  two-element  confocal  transducer. 
However,  in  response  to  the  reflected  beam  from  the  object,  the  transducer  acts  as  a  planar  reflector 
to  bounce  the  reflected  beam  back  toward  the  object.  This  field  (which  is  now  reflected  once  by 
the  object  and  once  by  the  transducer)  is  divergent.  Hence,  multiple  reflections  rapidly  loose  their 
strength  as  they  reflect  back  and  forth  between  the  object  and  the  transducer. 

2.a.4)  New  Accelerometer 

When  imaging  small  objects  (such  as  a  glass  bead  or  microcalcification)  in  a  viscoelastic 
background,  the  acoustic  emission  amplitude  can  be  small.  Detection  of  this  field  by  a  hydrophone 
would  not  be  too  difficult  if  the  frequency  is  high.  At  low  frequencies  however,  hydrophone  sensitivity 
usually  drops  to  a  level  that  may  not  be  sufficient  for  detection  of  such  small  particles.  As  an 
alternative  detector,  we  have  used  an  accelerometer.  The  new  accelerometer  (B&K,  Model  4387) 
has  a  useful  frequency  range  from  0.1  Hz  to  8  kHz,  with  a  sensitivity  of  26  mv/ms'2.  It  can  be 
placed  directly  in  contact  with  the  object  to  detect  small  vibrations.  A  potential  application  for  the 
accelerometer  is  detection  of  the  resonance  frequency  of  the  particles. 
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2.b  Phantom  Design 

To  test  the  performance  of  the  system  we  designed  three  phantoms. 

Phantom  1  is  made  from  RTV  (Room  Temperature  Vulcanizing)  material  and  includes  an 
embedded  glass  bead.  RTV  is  a  rubber-like  material  with  a  density  and  sound  velocity  close  to 
those  of  soft  tissues.  It  is,  however,  stiffer  than  most  soft  tissues.  The  attenuation  coefficient  of  this 
material  is  higher  than  most  soft  tissues.  The  glass  bead  with  a  diameter  of  400  pm  is  placed  at 
a  depth  of  3  mm  from  the  surface.  The  glass  bead  has  a  high  acoustic  impedance  similar  to  that 
of  breast  microcalcifications. 

Phantom  2  is  also  made  from  RTV,  but  it  includes  two  2-mm-diameter  holes  instead  of  the  glass 
bead.  The  holes  in  this  phantom  are  partially  air  filled,  representing  an  object  with  very  low  acoustic 
impedance.  Locations  of  the  holes  are  3  mm  and  10  mm  from  the  surface. 

Phantom  3  is  made  from  soft  gel  with  some  embedded  glass  beads  at  a  depth  of  10  mm.  The 
density  and  the  sound  speed  of  the  gel  material  is  very  close  to  that  of  soft  tissue,  and  its  stiffness 
is  similar  or  lower  than  many  soft  tissues.  This  material  has  an  attenuation  coefficient  much  lower 
than  that  of  soft  tissue. 

2.c  Acoustic  and  Vibrational  Noise  Reduction  in  the  Experimental  System 

2.C.1)  Quiet  Room 

The  quiet  room  is  built  within  the  Ultrasound  Research  Laboratory,  with  no  window  to  the 
outside.  Its  area  is  3x3  meters.  To  reduce  transmission  of  the  noise  from  other  rooms,  the  wall  and 
the  ceiling  is  covered  by  2  in  thick  Sonex  acoustical  foam.  The  foam  has  a  corrugated  wavy  surface 
reducing  reflection  of  sound.  The  room  hosts  the  air  table  and  the  scanning  tank  with  a  minimum 
of  electronic  equipment.  The  computer  and  other  electronic  instruments  that  have  fans  and  produce 
a  high  level  of  audio  noise  are  placed  outside  the  quiet  room  while  connected  via  cables  that  pass 
through  the  wall.  The  ambient  airborne  audio  range  noise  in  the  quiet  room  measures  approximately 
20  dB  less  than  the  ambient  noise  level  in  the  adjacent  laboratory. 

2.C.2)  Air  Table 

To  reduce  transmission  of  structural  vibrations  to  the  experiment  setup,  an  air  table  (Technical 
Manufacturing  Corp.,  Model  68-571)  is  now  used  to  support  the  water  tank  in  which  the  experiments 
take  place.  This  air  table  provides  vertical  isolation  efficiency  of  90-99%  at  10  Hz. 

2x3)  Reducing  the  Multipath  and  Reverberation  Effects —  New  Tank  Design  and  Absorbers 

Major  problems  that  led  us  to  design  a  new  water  tank  for  our  experiments  are  described  below 
(also  briefly  mentioned  in  Table  1): 

Multipath  problem  —  The  acoustic  emission  field  produced  as  a  result  of  object  vibration  in  a 
conventional  water  tank  can  travel  along  many  paths  before  reaching  the  hydrophone.  These 
field  components  can  interact  constructively  or  destructively  depending  on  the  path  lengths.  This 
phenomenon  causes  an  error  in  the  resulting  image  because  the  amplitude  of  the  recorded  acoustic 
emission  changes  with  respect  to  the  position  of  the  object  in  the  tank. 
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Reverberation  problem  —  When  the  acoustic  emission  field  travels  back  and  forth  between  two 
parallel  reflective  surfaces  of  the  tank,  the  acoustic  field  may  reach  a  resonance  state,  in  which 
the  field  amplitude  would  become  highly  dependent  on  the  wavelength  and  the  tank  size.  This 
phenomenon  can  particularly  be  a  problem  when  the  sound  reverberates  between  the  water  surface 
and  the  bottom  of  the  tank,  causing  the  acoustic  amplitude  to  vary  with  the  slightest  ripple  on  the 
water  surface. 

New  tank  design  —  To  reduce  both  the  multipath  and  reverberation  problems  a  new  tank  was 
designed  and  constructed.  The  new  tank  has  a  double  layer  bottom  with  the  top  layer  positioned  in 
a  15  degree  angle  with  respect  to  the  lower  layer.  The  interior  of  the  tank  is  covered  with  sound 
absorbing  material  (Sorbothane,  Sorbothane  Inc.).  This  material  is  highly  absorbent,  reducing  the 
multipath  and  reverberation  effects  in  the  tank.  The  angled  bottom  also  reduces  the  possibility  of 
reverberation  between  the  water  surface  and  the  tank  bottom. 

2.d  Improving  the  Scanning  and  Data  Acquisition  Hardware  and  Software  to  Increase  the 
Scanning  Speed  and  Decreasing  the  Acoustic  Noise  of  the  Scanning  Mechanism 

Slow  scanning  speed  problem  —  Because  our  imaging  method  operates  on  one  point  of  the 
object  at  a  time,  the  process  can  take  a  long  time.  In  addition,  due  to  the  mechanical  vibration  of 
the  scanning  mechanism,  the  scanner  needs  a  period  of  time  (usually  a  fraction  of  a  second)  to  settle 
after  it  moves  to  a  new  point.  This  time  plus  the  time  required  to  acquire  and  transfer  the  data  to 
the  disk  can  be  substantial. 

2.d.l)  High-speed  Scanning 

In  order  to  speed  up  the  scanning  process,  we  must  prevent  the  noise  and  vibrations  caused  by 
the  scanning  structure  from  interfering  with  the  acoustic  emission  data.  One  way  to  do  this  is  by 
increasing  the  acoustic  emission  frequency  beyond  the  spectral  bandwidth  of  the  scanner  generated 
noise,  such  that  the  noise  can  be  filtered  out  without  losing  the  acoustic  emission  signal. 

Based  on  this  strategy,  we  have  performed  high-speed  scanning  at  30-60  kHz.  A  programmable 
filter  was  used  to  remove  the  scanner  noise.  The  noise  bandwidth  in  our  system  is  about  20  kHz. 
In  this  method,  the  transmitting  transducer  continuously  transmits  the  ultrasound  while  the  scanner 
scans  the  object  without  stopping,  and  at  the  same  time  the  hydrophone  continuously  receives  the 
resulting  acoustic  emission  signal.  Using  this  method  an  object  can  be  scanned  in  a  fraction  of  the 
time  that  would  be  required  for  scanning  by  the  conventional  point-by-point  scanning.  The  resulting 
images  have  high  signal-to-noise  ratio.  The  only  drawback  of  this  method  is  that  it  can  not  be  used 
at  acoustic  emission  frequencies  lower  than  20  kHz. 

2.d.2)  Quiet-scanning  System 

The  step  motors  used  in  the  present  scanner  produce  loud  audio  noise  that  interferes  with  the 
acoustic  emission  produced  by  the  object.  To  reduce  the  vibration  and  acoustic  noise  of  the  scanning, 
a  new  scanning  system  is  being  designed  that  is  based  on  analog  servo  motors  instead  of  the  present 
step  motors.  The  new  scanning  system  has  optically  encoded  axes  for  precise  positioning.  This 
system  is  currently  under  construction. 
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2.d.3)  Tone  Burst  Scanning 

Tone  burst  scanning  is  a  method  we  have  used  to  avoid  the  ultrasound  standing  wave  effect  and 
acoustic  emission  multipath  effect  in  the  water  tank  (Table  1). 

Ultrasound  standing  wave  effect  is  caused  by  multiple  reflection  of  the  sound  between  a  highly 
reflective  object  and  the  transducer.  Acoustic  emission  multipath  effect  is  caused  by  multiple 
reflection  of  the  acoustic  emission  by  the  walls  or  other  reflectors  within  the  tank.  Both  these 
phenomena  reduce  the  quality  of  the  final  image  and  can  introduce  false  patterns  on  the  image. 

In  this  method  we  excite  each  transducer  element  with  tone  bursts  of  continuous  waves.  The 
length  of  the  bursts  (Tg)  is  chosen  to  be  at  least  several  times  the  acoustic  emission  period,  2ir/ Alo, 
but  short  enough  to  avoid  any  overlapping  with  its  multiple  reflections  from  the  nearby  structures 
such  as  the  tank  walls.  Assuming  that  r  represent  the  time  delay  between  the  acoustic  emission 
signal  and  its  reflection  coming  from  the  closest  structure  (which  could  be  the  transducer,  wall,  water 
surface,  etc.),  we  can  write: 


27t/A oo  «  Tb  «  r.  (1) 

There  are  two  major  drawbacks  for  this  method.  The  first  problem  is  limitation  on  acoustic 
emission  frequency.  To  satisfy  the  above  equation  we  often  have  to  use  a  high  acoustic  emission 
frequency,  such  that: 


Alo  »  27 r/r.  (2) 

In  our  perimental  setup  the  acoustic  emission  frequency  needs  to  be  about  40  kHz  or  more. 

The  second  problem  is  production  of  unwanted  acoustic  noise.  Applying  a  sudden  burst  of 
energy  to  a  transducer  usually  produces  a  transient  vibration  of  the  transducer  case.  This  vibration 
introduces  a  wideband  noise  in  the  entire  experimental  setup.  To  reduce  this  effect,  we  can  use  a 
gradual  ascending  and  descending  tone  burst,  which  is  described  in  the  next  section. 

2.d.4)  Gating  Tone  Bursts 

A  “Soft  Gate”  was  constructed  using  a  raised  cosine  shaping  circuit  to  reduce  rise  time  and  fall 
time  of  tone  burst  edges.  This  approach  greatly  reduces  the  acoustic  noise  from  transducer  element 
shock,  thus  improving  the  detection  threshold  of  the  system. 

To  test  this  method,  we  scanned  the  RTV  glass-bead  phantom  with  and  without  the  soft  gate. 
Without  this  gate,  the  noise  from  the  transducer  shock  overwhelmed  the  weak  acoustic  emission  from 
the  bead.  The  image  quality  was  greatly  improved  and  the  bead  became  readily  detectable  after  using 
the  soft  gate.  Experiments  on  the  RTV  phantom  with  holes  also  showed  the  same  improvement. 

2.d.5)  Aluminum  Free  Tank 

When  imaging  small  objects  such  as  glass  beads  (or  microcalcifications),  the  signal  level  is  very 
low,  hence  any  source  of  acoustic  noise  should  be  studied  and  hopefully  eliminated. 

We  discovered  that  pieces  made  from  aluminum  can  cause  a  considerable  acoustic  noise  when 
the  piece  is  placed  in  water.  Clamps,  rods,  and  other  structures  used  for  holding  the  transducer  and 
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the  object  in  the  water  tank  are  usually  made  from  aluminum.  Chemical  interaction  of  water  and 
aluminum  produces  a  “flickering”  sound  in  the  tank.  This  sound,  which  is  of  transient  nature  but 
repeated  in  random  pattern,  is  a  wideband  noise  which  can  not  be  filtered  out  at  the  frequency  range 
that  is  of  our  interest  for  acoustic  emission.  Its  amplitude  is  significant  compared  to  the  acoustic 
emission  from  small  objects. 

Eliminating  all  aluminum  from  the  water  tank  significantly  reduced  background  noise. 


2.e  Modeling  and  Evaluation  of  the  Acoustic  Frequency  Responses,  Q(Aco)  and  H( Aw) 

The  aim  of  our  method  is  to  image  an  object  based  on  its  mechanical  characteristics.  This  is 
achieved  by  vibrating  the  object  by  applying  a  highly  localized  oscillating  force  to  each  point  of  the 
object.  The  localized  force  is  produced  by  modulating  the  intensity,  and  thereby  the  radiation  force, 
of  an  ultrasound  field  at  low  frequencies  (normally  in  kHz  range).  The  resulting  sound  emitted  by 
the  object  is  a  function  of  object  mechanical  characteristics  and  the  location  of  the  excitation  point. 
The  image  is  produced  by  mapping  the  amplitude  or  phase  of  this  sound,  which  is  detected  by  a 
sensitive  hydrophone,  versus  position.  Figure  1  of  reference  [1]  (page  52)  illustrates  this  method. 

In  this  section  we  describe  the  relationship  between  the  image  and  the  properties  of  the  object 
and  the  medium.  More  detail  analysis  is  presented  in  references  [1],  [2],  and  [3]. 

Here  we  use  the  terminology  which  is  used  in  reference  [1]  and  call  this  method  “Ultrasound- 
Stimulated  Vibro-Acoustography”  (USVA).  We  will  also  use  the  notations  used  in  that  reference.  In 
particular,  we  note  that  Q(Aco)  is  represented  by  and  H(  Au>)  is  represented  by 

Consider  a  plane  ultrasound  beam  interacting  with  a  planar  object  of  zero  thickness  and  arbitrary 
shape  and  boundary  impedance  that  scatters  and  absorbs.  The  radiation  force  vector,  F,  arising  from 
this  interaction  has  a  component  in  the  beam  direction  and  another  transverse  to  it.  The  magnitude  of 
this  force  is  proportional  to  the  average  energy  density  of  the  incident  wa ve(E)  at  the  object,  where 
()  represents  the  time  average,  and  S,  the  area  of  the  projected  portion  of  the  object  [6] 

F  =  d  rS(E),  (3) 


where  dr  is  the  vector  drag  coefficient  with  a  component  in  the  incident  beam  direction  and  another 
transverse  to  it.  The  coefficient  dr  is  defined  per  unit  incident  energy  density  and  unit  projected 
area.  For  a  planar  object,  dr  is  numerically  equal  to  the  force  on  the  object.  Physically,  the  drag 
coefficient  represents  the  scattering  and  absorbing  properties  of  the  object  and  is  given  by  [6] 


dr 


+  ns 


/ 


7  coso:sdS 


7  sin  asdS , 


(4) 


where  p  and  q  are  the  unit  vectors  in  the  beam  direction  and  normal  to  it,  respectively.  The  quantities 
na  and  n,  are  the  total  absorbed  and  scattered  powers,  respectively,  and  7  is  the  scattered  intensity, 
all  expressed  per  unit  incident  intensity.  Also,  is  the  angle  between  the  incident  and  the  scattered 
intensity,  and  dS  is  the  area  element.  The  drag  coefficient  can  also  be  interpreted  as  the  ratio  of  the 
radiation  force  magnitude  on  a  given  object  to  the  corresponding  value  if  the  object  were  replaced 
by  a  totally  absorbing  object  of  similar  size.  This  is  because  |  dr  |  =  1  for  a  totally  absorbing 
object.  This  coefficient  can  be  determined  for  objects  of  different  shapes  and  sizes.  For  simplicity, 
we  assume  a  planar  object  normal  to  the  beam  axis.  In  this  case,  the  transverse  component  vanishes, 


10 


thus,  the  drag  coefficient  (force)  will  have  only  a  component  normal  to  the  target  surface  which  we 
denote  by  scalar  dr  ( F ). 

To  produce  a  dynamic  radiation  force,  one  can  use  an  amplitude  modulated  beam  [7].  Consider 
an  amplitude  modulated  incident  (ultrasonic)  pressure  field,  p(t),  as 

m  =  cos  (Awt/2)  cosuiQt,  (5) 


where  PUo,  Aw/2,  and  wo  are  the  pressure  amplitude,  modulating  frequency,  and  the  center  frequency, 
respectively.  In  our  analysis  and  experiments,  we  assume  that  the  condition  Aw  <<  wo  holds.  In 
such  a  case,  the  energy  density  of  the  incident  field  has  slow  variations  in  time.  To  discriminate 

the  slow  time  variations  of  a  function,  let  us  define  the  short-term  time  average  of  an  arbitrary 

t+Tj  2 

function  £(t)  over  the  interval  of  T  seconds  at  time  instance  t,  as  (£(/))T  =  y  f  £(r)dr,  which 


t-T/2 

is  a  function  of  t.  The  long-term  time  average  (or  simply,  the  time  average)  is  obtained  by  setting 
T  — ►  oo.  To  compute  the  short-term  time  average  of  the  acoustic  energy  density  relevant  to  field 
variations  at  Aw/2,  we  choose  T  longer  than  the  ultrasound  wave  period  but  much  shorter  than  the 
modulation  period,  that  is  27t/wo  <<  T  <<  47t/Aw.  Under  this  condition,  the  short-term  time 

average  ofp2(t)  is  (p2(t))T  =  %*-(l  +  cos  A cot).  The  energy  density  is  given  by  p2(t)/pc 2  where  p 
and  c  are  the  density  and  propagation  speed  in  the  medium  [8].  We  are  interested  in  the  time- varying 
component  of  the  short-term  time  average  of  the  energy  density.  Denoting  this  component  by  e^w(i), 

p2 

we  can  write:  e^o ,(/)  =  ^§z  cos  A  cot.  This  component  of  the  energy  density  produces  a  time  varying 
radiation  force  on  the  target  (Eq.  3)  at  frequency  Aw.  The  amplitude  of  this  force,  F&w,  is 


(6) 


This  equation  states  that  the  time-varying  force  amplitude  is  proportional  to  the  square  of  incident 
ultrasound  pressure,  or  equivalently,  to  the  incident  power.  If  the  object  moves  in  response  to  this 
force,  then  the  high-frequency  ultrasound  energy  would  convert  to  low-frequency  mechanical  energy. 

The  radiation  force  F&w  vibrates  the  target  object  at  frequency  Aw.  Object  vibration  results  in 
an  acoustic  field  in  the  medium  (acoustic  emission).  This  field  is  related  to  object  shape,  size,  and 
viscoelastic  properties.  To  present  a  conclusive  analysis  of  this  relationship,  we  have  to  assume  an 
object  with  specific  characteristics.  Here,  we  assume  that  the  vibrating  object  has  a  circular  cross 
section  of  radius  b  and  uniformly  vibrates  back  and  forth  like  a  piston.  This  choice  allows  us  to 
illustrate  the  concept  in  a  simple  form.  We  also  consider  an  area  S  <  i rb2  of  the  piston  surface  to 
be  projected  normally  by  the  beam.  Similar  solutions  can  be  carried  out  for  other  objects. 

The  steady  state  normal  velocity  amplitude  of  a  piston,  Uau,  due  to  a  harmonic  force  Fau  at 
frequency  Aw,  can  be  described  in  terms  of  the  mechanical  impedance  Zau>, 


U Auj 


Faqj 

w 


(7) 


where  Zaui  =  Z'm  +  Zr  is  comprised  of  the  mechanical  impedance  of  the  object  in  vacuum  Z'm,  and 
the  radiation  impedance  of  the  object  Zr,  all  defined  at  Aw.  Modeling  the  object  as  a  mass-spring 
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(t 


system,  Z'm  can  be  written  in  terms  of  Au ;  as 

Z'm  =  R'm-  3  .  (8) 

where  m,  R'm,  and  K'  are  the  mass,  mechanical  resistance,  and  the  spring  constants  of  the  object, 
respectively.  The  radiation  impedance  of  the  piston  can  be  written  as  [9] 


where 

Zr  =  7T  b2(Rr  —  jXr), 

(9) 

Awf>Jl(2AaJ]’ 

(10) 

and 

7 r 

4 pc  f  .  f2Aub  \  .  2  , 

Xr  = -  /  sin  - cos  ol  }  sm  aaa, 

K  J  \  C  J 

(11) 

o 


where  Ji(-)  is  the  first-order  Bessel  function  of  the  first  kind.  In  many  applications  of  our  interest, 
the  wavelength  is  much  greater  than  the  object  size,  hence,  ^rb  —>  0.  In  such  cases  Zr  assumes 
a  simpler  form: 


The  mechanical  impedance  of  the  piston  object  can  now  be  written  as 
Zaw  =  ( R'm  +  nb^Rr)  -  j  fmAoj  -  ~  +  wb2Xr'\ 

^  u  f  mi 

7T p _ 4  .  2 \  ./  A  K1  8  i3a  \  Aub  n  K  } 

+  ^  3^  )’ 

Once  we  calculate  Uauj,  we  can  calculate  the  pressure  field  it  produces  in  the  medium.  We  assume 
that  the  acoustic  emission  signal  propagates  in  a  free  and  homogenous  medium.  The  farfield  acoustic 
pressure  due  to  a  piston  source  of  radius  b  set  in  a  planar  boundary  of  infinite  extent  is  given  by  [9], 


Pau,  =  -jAojp 


exp  (jAuil/c) 
47 tI 


2Ji{ 


Autb 

c 


sin  7?) 


Aw  6 


sin  7? 


x 


COST? 


COST?  +  /3jg 


{2WUau), 


(14) 


where  l  is  the  distance  from  the  observation  point  to  the  center  of  the  piston,  t?  is  the  angle  between 
this  line  and  the  piston  axis,  and  (3b  is  the  specific  acoustic  admittance  of  the  boundary  surface. 
The  specific  acoustic  admittance  is  (3b  —  where  Zb,  the  acoustic  impedance  of  the  boundary, 
represents  the  ratio  between  the  pressure  and  normal  fluid  velocity  at  a  point  on  the  surface.  The 
factor  of  two  comes  from  the  presence  of  the  boundary  wall.  It  would  be  replaced  by  unity  if  the 
boundary  wall  were  not  present  [8].  The  acoustic  emission  field  resulting  from  object  vibration  can 
be  written  in  terms  of  the  incident  ultrasound  pressure  by  combining  Eqs.  6,  7,  and  14,  as 


P/Alj  =  {  j 


.  A  lu  exp  (j  A u)  l  /  c) 


47 d 


2Ji(— ^sinr?)  cost? 

X 


Aujb  {)  "  cos  3)  +  (3b 


sin' 


x 


[(R'm  +  TTb2Rr)  -  j  (mAu  -j£;  +  -Kb2Xr)\ 


(27 Tb2)P20Sdr, 


(15) 
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For  wavelengths  that  are  long  compared  to  the  object  size,  i.e.,  when  bAu/c  — ►  0,  the  term  in 
the  first  brace  approaches  a  constant,  hence  we  may  consider  the  contents  of  the  first  brace  to  be  an 
object  independent  function  (the  specific  acoustic  admittance  relates  to  the  surrounding  boundary 
surface).  Under  these  conditions,  the  first  brace  in  the  above  equation  represents  the  effect  of  the 
medium  on  the  acoustic  emission  field,  which  we  may  call  the  medium  transfer  function,  and  denote 
it  by 

2J1(^sintf)  cos'd 

+  (16) 

The  second  bracket  in  Eq.  15  is  or  the  mechanical  admittance  of  the  object  at  the  frequency  of 
the  acoustic  emission  (Aw),  and  we  denote  it  by  Y&w.  It  is  convenient  to  combine  this  term  with  the 
next  term  ( 2vb 2)  in  Eq.  15,  as  Qau>  =  27rf>2YAo>  =  2?r if  / Zaw>  which  is  the  total  acoustic  outflow1  by 
the  object  per  unit  force.  Function  Qau  represents  the  object  characteristics  at  the  acoustic  frequency. 
We  may  thus  rewrite  Eq.  15  in  a  more  compact  form  as 

Paw  =  HAu(l)QAoJP%oSdr.  (17) 

Equation  17  indicates  that  the  acoustic  emission  pressure  is  proportional  to:  1)  the  square  of  ultrasound 
pressure  Puo,  2)  the  ultrasound  characteristics  of  the  object,  dr,  in  the  projected  area  S,  3)  the 
acoustic  outflow  by  this  object,  Qau,  representing  the  object  size  b  and  its  mechanical  admittance  at 
the  acoustic  frequency,  VAa),  and  4)  the  transfer  function  of  the  medium  at  the  acoustic  frequency. 
The  above  equation  illustrates  the  basic  nonlinear  relationship  between  the  ultrasound  and 
acoustic  emission  pressure  amplitudes.  Note  that  neither  the  medium  nor  the  object  need  to  be 
nonlinear  for  this  relationship  to  hold.  It  is  interesting  to  note  that  the  projection  area  S  and  the 

vibrating  area  Tib2  play  different  roles.  The  projection  area  determines  the  extent  of  the  force  applied 

to  the  object  (Eq.  6).  The  vibrating  area,  however,  influences  the  total  acoustic  outflow  in  the  medium 
caused  by  object  vibration.  The  mechanism  of  object  vibration  is  somewhat  analogous  to  that  of  a 
loudspeaker,  where  the  electromotive  force  is  exerted  to  a  small  area  of  the  membrane  (usually  at 
the  center)  causing  the  entire  membrane  surface  to  vibrate.  In  our  method,  the  size  of  the  vibrating 
area  depends  on  the  object  structure.  For  a  free  suspended  point  object,  smaller  than  the  beam  cross 
section,  the  vibrating  area  would  be  the  same  as  the  projection  area.  For  a  large  stiff  plate,  however, 
the  vibrating  area  could  be  much  larger  than  the  projected  area  (similar  to  a  loudspeaker).  In  some 
cases  it  is  more  convenient  to  have  the  acoustic  emission  field  in  terms  of  the  applied  force  F^. 
Referring  to  Eq.  6,  we  can  rewrite  Eq.  17  as 


CL  47 ll 


Paiv  4pc  H Au;(l)Q Auj P Atj •  (18) 

Again  in  analogy  to  a  loudspeaker,  i'Aw,QAw>and  #Aw(0  represent  the  electromotive  force, 
dynamic  characteristics  of  the  membrane,  and  the  propagation  medium  transfer  function. 


Acoustic  outflow  is  the  volume  of  the  medium  (e.g.,  the  fluid)  in  front  of  the  object  surface  that  is  displaced  per  uni;  time  due  to  object  vibration. 
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2.f  Beam  Forming:  Modeling  and  Experimental  Results 


Beam  forming  for  our  method  is  basically  different  from  that  used  in  the  conventional  ultrasound 
system.  In  a  conventional  pulse-echo  system  beam  forming  is  done  to  produce  the  smallest  possible 
beam  width  with  the  highest  intensity.  Here  we  are  interested  in  producing  a  radiation  stress  field 
in  a  small  region  in  a  three-dimensional  space. 

We  have  developed  three  beam  forming  methods:  (1)  single  modulated  beam,  (2)  two  crossing 
beam  (x-focal),  and  (3)  the  confocal  beams.  The  single  beam  is  of  theoretical  interest  and  is  not 
used  in  our  practice  because  of  its  low  depth  selectivity.  The  theory  of  confocal  beam  is  explained 
in  reference  [1],  The  theory  and  experimental  results  for  the  x-focal  beam  and  a  comparison  with 
the  other  two  beams  is  presented  in  [4].  Here  we  mainly  present  the  theory  of  the  confocal  beam, 
which  is  of  the  most  interest  in  our  experimental  work. 

Confocal  beam  is  produced  by  projecting  two  coaxial,  confocal,  continuous  wave  (CW)  ultra¬ 
sound  beams  on  the  object.  An  amplitude  modulated  field  is  produced  only  at  the  interference  region 
of  the  two  unmodulated  beams  around  their  focal  areas,  resulting  in  a  small  resolution  cell.  For  this 
purpose,  elements  of  a  two-element  spherically  focused  annular  array  (consisting  of  a  central  disc 
with  radius  ai  and  an  outer  ring  with  the  inner  radius  oia!2  and  outer  radius  of  <22)  are  excited  by 
separate  CW  signals  at  frequencies  =  u>q  -  Au>/2  and  022  =  ujq+Auj/2.  We  assume  that  the  beams 
are  propagating  in  a  lossless  medium,  in  the  +z  direction  of  a  Cartesian  coordinate  system  (x,y,  z), 
with  the  joint  focal  point  at  2  =  0.  The  resultant  pressure  field  on  the  z  =  0  plane  may  be  written  as 

P(t)  =  Pi(r)cos(uJit  +  il>i(r))  +  P2(r)cos(u}2t  +  i>2{r)),  (19) 


where  r  =  yjx2  +  y2  is  the  radial  distance.  The  amplitude  functions  are  [8,  10] 


7r  at 


Pi(r)  =  pcUoi- — -  jinc 


Aizo 


and 


'  rai  \ 

Mz0J' 


P2{r)  =  pcUo2 


7T 


^2Z0 


rjmc(S)-aMS)j' 


(20) 


(21) 


7t~r 
Zq  ' 


where  A t-  =  27r/u;t-,  i=l,  2,  is  the  ultrasound  wavelength,  Uoi  is  the  particle  velocity  amplitude  at  the 
i-th  transducer  element  surface,  and  jinc  (X)  =  Ji(2-kX)/-kX.  The  phase  functions,  ^(r)  =  -fj 
for  i=l,  2,  are  conveniently  set  to  be  zero  at  the  origin. 

Now,  we  define  a  unit  point  target  at  position  (m0, 2/0)  on  the  focal  plane  with  a  drag  coefficient 
distribution  as 


dr{x,  y)  =  6(x  -x0,y-  y0 ),  (22) 

such  that  dr(x,  y)dxdy  is  unity  at  the  (x0,y0)  and  zero  elsewhere.  This  equation  is  merely  used  as  a 
mathematical  model  because  dT  is  physically  finite.  In  this  case,  the  projected  area  can  be  considered 
to  be  S  =  dxdy.  We  replace  drS  in  Eq.  3  with  dr(x,y)dxdy,  and  follow  the  steps  similar  to  those 
outlined  in  Eqs.  5  and  6  for  the  pressure  field  expressed  by  Eq.  19,  then  the  complex  amplitude  of 
the  normal  component  of  the  force  on  the  unit  point  target  can  be  found  as 

roai\  [  7ra|  ( rpt^  _  ?ra£  ..n  /rpOg^ 


FaU^Vo)  =  \pU0iUo2j^ 

r%Au) 


-Jinc 


( 


A12&/ 


x  exp 


\  2cz0  ) 


(23) 
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where  the  arguments  xo&ndyo  are  added  to  denote  the  position  of  the  point  target,  and  ro  = 
\jx o  +  Vq-  Equation  23  describes  the  spatial  distribution  of  the  force  (the  stress  field).  This  equation 
shows  that  the  stress  field  is  confined  to  the  regions  near  the  beam  axis  (ro  =  0),  and  decays  as 
the  radial  distance  tq  increases.  The  lateral  extent  of  the  stress  field,  and  hence  the  resolution  cell 
diameter,  would  be  smaller  at  higher  ultrasound  frequencies  (smaller  Ai  and  A2).  One  can  calculate 
the  total  force  on  an  arbitrary  object  by  integrating  the  force  over  the  projected  area.  The  axial  extent 
of  the  resolution  cell  (depth  resolution,  or  the  depth  of  field)  can  be  determined  by  calculating  the 
force  jFaw  as  a  function  of  the  depth  variable,  in  a  fashion  as  outlined  in  Eqs.  19  to  23. 

2.g  Point  Spread  Function 

To  produce  an  image  we  scan  the  object  in  a  plane  and  record  the  complex  amplitude  of  the 
acoustic  emission,  Paw.  at  different  positions.  In  this  process,  we  keep  Auj  fixed.  For  transverse 
view  images,  the  scan  plane  is  the  focal  plane  (x  —  y).  Alternatively,  for  the  parallel  view  the  scan 
plane  is  the  x  —  z  plane.  In  the  conventional  ultrasound  imaging  context,  these  two  views  are  called 
the  C-scan  and  B-scan,  respectively.  Our  main  focus  here  is  the  transverse  view  imaging.  In  this 
case,  the  acoustic  emission  data  obtained  by  vibrating  the  object  at  point  (x,y)  are  assigned  to  the 
corresponding  point  (x,y)  in  the  image. 

Before  defining  the  image,  we  need  to  define  the  function  that  represents  the  object.  Referring 
to  Eq.  17,  the  terms  that  are  object  dependent  are  the  drag  coefficient  dT  and  the  function  Qaw 
(assuming  that  H^^l)  is  object  independent).  The  object  function  g(x,y)  is  thus  defined  as  the 
spatial  distribution  of  these  terms, 


g(x,y)  =  Q&u>(x,y)dr{x,y).  (24) 

Variables  x  and  y  are  added  to  denote  the  dependency  of  dr  and  Qaw  on  position.  In  particular, 
QAw(x,y)  implies  the  total  acoustic  outflow  by  the  object  when  the  force  is  applied  at  point  (x,y). 

Commonly,  an  imaging  system  is  studied  through  its  point  spread  function  (PSF),  which  is 
defined  as  the  image  of  a  point  object.  To  determine  the  PSF  of  our  system,  we  consider  a  unit  point 
target  at  the  origin  with  unit  mechanical  response,  Q&.w(x,y)  =  1.  Hence,  referring  to  Eqs.  22  and 
24,  we  can  write  g(x,y)  =  6(x,y).  To  obtain  the  PSF,  we  move  this  point  object  to  every  possible 
position  (xo,i/o)  on  the  z  =  0  plane  and  form  the  image  using  the  resulting  acoustic  emission  field, 
Pauj(xo,  Vo)-  Since  so  and  yo  are  now  being  treated  as  variables,  we  may  replace  them  by  variables 
x  and  y,  respectively.  We  define  the  normalized  PSF  of  the  coherent  imaging  system  as  the  complex 
function 

h(x,y)  =  PAu(x,y)/P^(0,0).  (25) 

Division  by  Pao>(0,  0)  cancels  the  constant  multipliers.  Referring  to  Eqs.  18  and  23,  we  can  write 


This  equation  illustrates  that  the  system  PSF  is  a  circularly  symmetric  function  with  the  peak  at  the 
origin  and  decaying  amplitude  with  increasing  radial  distance  r.  Amplitude  decays  faster  for  higher 
ultrasound  frequency.  This  function  will  be  discussed  further  in  the  next  section. 
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where  the  arguments  x§  and  yo  are  added  to  denote  the  position  of  the  point  target,  and  ro  = 
\Jx q  +  Vq.  Equation  23  describes  the  spatial  distribution  of  the  force  (the  stress  field).  This  equation 
shows  that  the  stress  field  is  confined  to  the  regions  near  the  beam  axis  (ro  =  0),  and  decays  as 
the  radial  distance  ro  increases.  The  lateral  extent  of  the  stress  field,  and  hence  the  resolution  cell 
diameter,  would  be  smaller  at  higher  ultrasound  frequencies  (smaller  Ai  and  A2).  One  can  calculate 
the  total  force  on  an  arbitrary  object  by  integrating  the  force  over  the  projected  area.  The  axial  extent 
of  the  resolution  cell  (depth  resolution,  or  the  depth  of  field)  can  be  determined  by  calculating  the 
force  F&w  as  a  function  of  the  depth  variable,  in  a  fashion  as  outlined  in  Eqs.  19  to  23. 

2.g  Point  Spread  Function 

To  produce  an  image  we  scan  the  object  in  a  plane  and  record  the  complex  amplitude  of  the 
acoustic  emission,  P^,  at  different  positions.  In  this  process,  we  keep  Au;  fixed.  For  transverse 
view  images,  the  scan  plane  is  the  focal  plane  (x  —  y).  Alternatively,  for  the  parallel  view  the  scan 
plane  is  th  ex  —  z  plane.  In  the  conventional  ultrasound  imaging  context,  these  two  views  are  called 
the  C-scan  and  B-scan,  respectively.  Our  main  focus  here  is  the  transverse  view  imaging.  In  this 
case,  the  acoustic  emission  data  obtained  by  vibrating  the  object  at  point  (x,y)  are  assigned  to  the 
corresponding  point  (x,y)  in  the  image. 

Before  defining  the  image,  we  need  to  define  the  function  that  represents  the  object.  Referring 
to  Eq.  17,  the  terms  that  are  object  dependent  are  the  drag  coefficient  dr  and  the  function  Q&w 
(assuming  that  H&w(l)  is  object  independent).  The  object  function  g(x,y)  is  thus  defined  as  the 
spatial  distribution  of  these  terms, 


g{x,  y)  =  Q&u{x,y)dr(x,y).  (24) 

Variables  x  and  y  are  added  to  denote  the  dependency  of  dr  and  Q/\,,,  on  position.  In  particular, 
Qaw(x,v)  implies  the  total  acoustic  outflow  by  the  object  when  the  force  is  applied  at  point  (x,y). 

Commonly,  an  imaging  system  is  studied  through  its  point  spread  function  (PSF),  which  is 
defined  as  the  image  of  a  point  object.  To  determine  the  PSF  of  our  system,  we  consider  a  unit  point 
target  at  the  origin  with  unit  mechanical  response,  Qaoj(x-,  y)  =  1.  Hence,  referring  to  Eqs.  22  and 
24,  we  can  write  g(x,y )  =  S(x,y).  To  obtain  the  PSF,  we  move  this  point  object  to  every  possible 
position  ( xo,yo )  on  the  2:  =  0  plane  and  form  the  image  using  the  resulting  acoustic  emission  field, 
Pau>{x0i  yo)-  Since  xq  and  yo  are  now  being  treated  as  variables,  we  may  replace  them  by  variables 
x  and  y,  respectively.  We  define  the  normalized  PSF  of  the  coherent  imaging  system  as  the  complex 
function 

h{x,y)  =  PAu,(x,y)/P&u(0,0).  (25) 

Division  by  Paw(0,  0)  cancels  the  constant  multipliers.  Referring  to  Eqs.  18  and  23,  we  can  write 


This  equation  illustrates  that  the  system  PSF  is  a  circularly  symmetric  function  with  the  peak  at  the 
origin  and  decaying  amplitude  with  increasing  radial  distance  r.  Amplitude  decays  faster  for  higher 
ultrasound  frequency.  This  function  will  be  discussed  further  in  the  next  section. 
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2.j  Discussion 

The  theory  presented  here  clarifies  the  fundamental  relationships  between  the  object  parameters 
and  its  USAE  image.  This  theory  indicates  that  the  acoustic  emission  from  an  object  is  related  to 
its  mechanical  admittance  and  the  size  of  the  object  (which  acts  as  a  source  of  acoustic  emission). 
Noting  that  microcalcifications  are  very  small  in  size,  detecting  the  acoustic  emission  field  from  such 
a  small  source  becomes  a  challenge.  Hence,  to  improve  the  detectability,  the  acoustic  noise  in  the 
medium  must  be  minimized.  This  includes  the  ambient  noise,  noise  from  the  scanner,  vibrations 
produced  by  the  transducer  body,  etc.  The  improvements  that  we  have  made  in  the  laboratory  setup, 
including  using  the  quiet  room,  air  table,  absorbing  material,  soft  gate,  etc.,  have  enabled  us  to 
considerably  increase  the  signal-to-noise  ratio  of  our  data.  Now,  with  our  low  noise  setup,  we  can 
detect  glass  beads  in  RTV  material.  This  was  not  possible  before  making  these  improvements.  These 
results  are  in  support  of  the  initial  hypothesis  that  USAE  can  be  used  to  detect  particles  of  small  size 
in  a  viscoelastic  material.  Using  the  tone  burst  method  has  enabled  us  to  eliminate  the  standing  wave 
problem  AT  the  expense  of  working  at  a  higher  acoustic  emission  frequency.  We  need  to  investigate 
other  methods,  such  as  using  a  flat  annular  array  transducer,  to  be  able  to  reduce  ultrasound  standing 
waves  when  working  at  lower  acoustic  emission  frequencies.  Driving  a  particle  at  its  resonance 
could  enhance  detectability  and  selectivity.  That  is  because  at  resonance,  particle  motion  is  at  its 
maximum.  Also,  because  the  resonance  frequency  of  a  particle  depends  on  its  size  and  mass,  this 
method  can  be  used  to  selectively  detect  particles  of  particular  size  and  mass.  Particle  resonance  is 
another  subject  of  interest  that  remains  to  be  studied  further. 

2.k  Key  Research  Accomplishments 

•  Development  of  the  theoretical  model  describing  the  relation  between  the  acoustic  emission  field 
and  the  mechanical  parameters  of  the  object,  as  well  as  to  the  parameters  of  the  propagation 
media. 

•  Developing  the  theory  describing  the  point  spread  function  of  the  imaging  system  with  three 
different  beamforming  methods. 

•  Improvement  of  the  sensitivity  of  the  system  to  detect  small  particles  in  a  viscoelastic  medium. 
Now  we  can  detect  a  glass  bead  with  diameter  of  260  /xm  (mimicking  a  microcalcification)  in 
gel  (mimicking  soft  tissue). 

•  Identifying  causes  of  image  artifacts,  such  as  ultrasound  standing  wave,  multipath,  and  rever¬ 
beration  of  acoustic  emission. 

•  Developing  methods  for  reducing  or  eliminating  acoustic  emission  multipath  and  reverberation 
problems. 

•  Developing  methods  for  reducing  acoustic  and  vibrational  noise  in  the  system 

2.1  Reportable  Outcomes 

1.  Fatemi,  M.,  and  J.  F.  Greenleaf:  Ultrasound-stimulated  vibro-acoustic  spectrography.  Science 
280:82-85,  April  3,  1998  (pages  24-27  of  this  report). 

2.  Fatemi,  M.,  and  J.  F.  Greenleaf:  Vibro-acoustography:  An  imaging  modality  based  on  ultrasound- 
stimulated  acoustic  emission,  proceedings  of  the  National  Academy  of  Sciences  USA  (In  Press) 
(pages  28-56  of  this  report). 
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3.  Fatemi,  M.,  and  J.  F.  Greenleaf:  Potential  applications  of  ultrasound-stimulated  vibro- 
acoustography  in  medical  diagnosis  and  material  characterization.  137th  Meeting  of  the  Acous¬ 
tical  Society  of  America  and  2nd  Convention  of  the  European  Acoustical  Association,  Berlin, 
Germany,  March  14-19,  1999  (Invited  presentation)  (page  58  of  this  report). 

4.  Fatemi,  M.,  and  J.  F.  Greenleaf:  Beam  forming  for  ultrasound-stimulated  vibro-acoustography. 
ICTCA’99  Fourth  International  Conference  on  Theoretical  and  Computational  Acoustics,  Trieste, 
Italy,  May  10-14,  1999  (Keynote  presentation)  (page  59  of  this  report). 

2.m  Conclusions 

In  this  report  we  presented  analytical  relationships  that  relate  the  acoustic  emission  field  to  the 
mechanical  parameters  of  object.  This  model  provides  a  tool  for  better  understanding  of  the  USAE 
imaging,  as  well  as  for  characterizing  objects  in  terms  of  their  mechanical  parameters.  Noting  that 
the  mechanical  stiffness  of  breast  microcalcification  is  drastically  different  from  that  of  soft  tissue, 
characterizing  mirocalcifications  based  on  their  mechanical  parameters  is  of  particular  interest.  We 
have  also  presented  the  point  spread  function  of  the  imaging  system  by  which  we  can  model  the 
system.  We  have  shown  that  the  experimental  system  can  be  improved  to  be  able  to  image  small 
particles,  about  the  size  of  breast  microcalcifications,  in  a  viscoelastic  medium. 

2.m.l)  So  what? 

We  have  shown  that  ultrasound-stimulated  acoustic  emission  method  (or  vibro-acoustography)  is 
capable  of  detecting  small  particles  in  a  viscoelastic  medium.  Improvement  of  this  technique  such 
that  it  can  be  used  to  reliably  detect  breast  microcalcifications  would  be  a  significant  contribution 
for  breast  cancer  detection.  This  tool  may  particularly  be  of  interest  in  examining  dense  breasts,  or 
pregnant  patients,  where  x-ray  mammography  cannot  be  used. 
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Figures  and  Legends 


APPENDIX  A 


Figure  1  Image  of  glass  beads  in  gel.  This  image  is  obtained  by  scanning  the  gel  phantom.  This  image  covers  an  area  of  20  mm 
by  20  mm,  scanned  at  0.2  mm/pixel.  The  beads  are  located  about  10  mm  deep  inside  the  gel.  The  diameter  of  the  top  bead  is  about 
260  ft. m,  and  the  diameter  of  the  largest  bead  is  about  400  fx m.  The  image  is  obtained  by  tone  burst  method  and  using  the  soft  gate. 
The  acoustic  emission  frequency  was  set  at  53.5  KHz. 


Figure  2  Image  of  a  glass  bead  in  RTV.  This  image  is  obtained  by  scanning  the  RTV  phantom  containing  glass  bead.  The  image 
covers  an  area  of  10x10  mm  and  is  scanned  at  0.1  mm  increments.  The  glass  bead  is  located  at  3  mm  from  the  surface  of  the 
phantom,  and  its  diameter  is  about  350  fim.  The  image  is  obtained  by  tone  burst  method  and  using  the  soft  gate.  The  acoustic 
emission  frequency  was  set  at  46.6  KHz. 
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Figure  3  Image  of  holes  in  RTV.  This  image  covers  a  square  of  60x60  mm,  and  is  scanned  at  0.3  mm/pixel.  The  two  holes  (diameters 
about  2  mm)  are  located  at  3  mm  and  10  mm  from  the  surface.  The  ultrasound  beam  was  focused  near  the  surface  of  the  phantom. 
The  image  is  obtained  by  tone  burst  method  and  using  the  soft  gate.  The  acoustic  emission  frequency  was  set  at  50  KHz. 
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Ultrasound-Stimulated  Vibro-Acoustic 
Spectrography 

Mostafa  Fatemi  and  James  F.  Greenleaf 

An  ultrasound  method  based  on  radiation  force  is  presented  for  imaging  the  acoustic 
response  of  a  material  to  mechanical  excitation.  Acoustic  energy  was  emitted  from  solids 
and  tissues  in  response  to  an  oscillatory  radiation  force  produced  by  interfering  focused 
beams  of  ultrasound.  Frequency  spectra  of  ultrasound-stimulated  acoustic  emission 
exhibited  object  resonances.  Raster-scanning  the  radiation  force  over  the  object  and 
recording  the  amplitude  and  phase  of  the  emitted  sound  resulted  in  data  from  which 
images  related  to  the  elastic  compositions  of  the  acoustically  emitting  objects  could  be 
computed.  Acoustic  emission  signals  distinguished  tuning-fork  resonances,  submilli¬ 
meter  glass  spheres,  and  calcification  in  excised  arteries  and  detected  object  motions 
on  the  order  of  nanometers. 


The  mechanical  response  of  objects  to  ex¬ 
ternal  forces  is  of  considerable  interest  in 
medical  diagnosis,  nondestructive  inspec¬ 
tion  of  materials,  and  materials  science.  An 
applied  force  is  often  used  to  produce  dis¬ 
placement  from  which  elastic  constants, 
like  spring  constants,  can  be  determined.  In 
resonant  ultrasound  spectroscopy,  an  ultra¬ 
sound  source  and  a  detector  are  used  to 
measure  the  resonance  frequencies  of  a  sam¬ 
ple  with  known  size  and  mass.  The  reso¬ 
nances  are  related  to  mechanical  parame¬ 
ters,  including  the  elastic  constants  of  the 
material  (1).  Recently,  a  magnetic  reso¬ 
nance  elastography  method  for  quantita¬ 
tively  measuring  the  displacement  of  tissues 
in  response  to  externally  applied  cyclic  forc¬ 
es  was  reported  by  Muthupillai  et  al.  (2). 
The  method  resulted  in  high-resolution  im¬ 
ages  of  the  shear  modulus  of  normal  and 
pathologic  tissues.  Others  have  used  ultra¬ 
sound  to  measure  tissue  displacement  asso¬ 
ciated  with  externally  applied  compressive 
and  cyclic  forces  (3). 

We  describe  an  imaging  technique  that 
uses  acoustic  emission  to  map  the  mechanical 
response  of  an  object  to  local  cyclic  radiation 
forces  produced  by  interfering  ultrasound 
beams.  Radiation  force  is  generated  by 
changes  in  the  energy  density  of  an  acoustic 
field  (4).  For  instance,  a  collimated  ultra¬ 
sound  beam  impinging  normally  on  the  sur¬ 
face  of  an  object  of  arbitrary  shape  and 
boundary  impedance  will  produce  a  radiation 
force.  The  radiation  force  arising  from  this 
interaction  has  a  component  F  =  drs( E)  (5) 
in  the  beam  direction.  This  component  is 
proportional  to  the  time-average  energy  den¬ 
sity  of  the  incident  wave  (E),  the  projected 
area  of  the  object  s,  and  dr  (6),  the  scattering 
and  absorbing  properties  of  the  object. 

We  probe  the  object  by  arranging  the 
intersection  of  two  focused  continuous- 


wave  (CW)  ultrasound  beams  of  different 
frequencies  at  a  selected  point  on  the  ob¬ 
ject.  Interference  in  the  intersection  region 
of  the  two  beams  produces  sinusoidal  mod¬ 
ulation  of  the  ultrasound  energy  density. 
Modulation  of  the  energy  density  creates  an 
oscillatory  force,  effectively  vibrating  the 
object  at  the  selected  region.  The  resulting 
vibration  of  the  object  produces  an  acoustic 
field  [acoustic  emission  (7)]  that  can  be 
measured  some  distance  away. 

Ultrasound  beams  can  be  constructed  in 
a  variety  of  ways  for  this  purpose  (8).  We 
used  two  coaxial,  confocal  transducer  ele¬ 
ments  of  a  spherically  focused  annular  array 
(consisting  of  a  central  disc  and  an  outer 
annulus)  driven  by  two  CW  signals  at 
slightly  different  frequencies  o>1  and  w2 
(Fig.  1).  The  energy  density  at  a  point  in 
this  ultrasound  field,  say  at  the  focus,  is 
proportional  to  the  square  of  the  sum  of  the 
ultrasound  fields  from  the  two  elements. 
Squaring  the  sum  of  two  sines  gives  rise  to 
sum  and  difference  frequency  terms.  Thus, 


high-frequency  and  low-frequency  varia¬ 
tions  in  energy  density  result  at  the  inter¬ 
section  of  the  two  beams  produced  by  the 
two  elements.  Ultrasound-stimulated  acous¬ 
tic  emission  results  from  the  energy  term 
that  produces  a  low-frequency  vibration. 
The  low-frequency  force  on  a  target  at  the 
focal  point  can  be  computed  by 

F,(t)  =  drJ  |  (EfxJt,x,y))Tdxdy 
s 

=  Cdr  cos(Awt)  (1) 

where  C  is  a  constant,  Aw  =  |  Wj  -  w2| ,  S 
is  the  area  over  which  Efocal(t,x,y),  the  total 
energy  density  in  the  focal  plane,  has  sig¬ 
nificant  value,  and  <)T  represents  a  short¬ 
term  time  average  (9).  For  focused  beams, 
the  intersection  region  can  be  small  enough 
that  Fj(t)  can  be  considered  to  be  an  oscil¬ 
lating  point  force  applied  to  the  object  at 
the  focal  intersection  of  the  beams. 

To  produce  an  ultrasound-stimulated  vi- 
bro-acoustic  spectrogram,  we  vibrate  a  small 
region  of  the  object  with  an  oscillating 
radiation  force  of  varying  frequency.  The 
complex  amplitude  of  the  resulting  acoustic 
emission  field  is 

*(Aw)  =  CdrH(Aw)Q(Aw)  (2) 

where  Q(Aw)  is  a  complex  function  repre¬ 
senting  the  mechanical  frequency  response, 
or  admittance,  of  the  object  at  the  selected 
point,  and  H(Aw)  represents  the  combined 
frequency  response,  or  transfer  function,  of 
the  propagation  medium  and  receiver  and  is 
assumed  to  be  fixed  and  known  (10).  Re¬ 
cording  <£(Aw)  allows  us  to  obtain  Q(Aw) 
for  each  point  within  a  constant  multiplier 
(II).  We  raster-scan  the  radiation  force 
over  the  object  to  produce  data,  which  can 


Fig.  1.  Experimental 
system  for  ultrasound- 
stimulated  vibro-acous- 
tic  spectrography:  a  two- 
element  confocal  ultra¬ 
sound  annular  array 
transducer,  consisting  of 
a  center  disc  and  an  out¬ 
er  ring.  The  elements  are 
driven  by  two  CW  sourc¬ 
es,  at  frequencies  equal 
to  co1  and  w2  =  w1  +  Aw, 
where  these  frequencies 

are  very  close  to  the  central  frequency  of  the  elements,  and  Aw  is  much  smaller  than  (<1  %)  the  center 
frequency  of  the  ultrasound  transducer.  The  beams  interact  only  in  a  small  region  around  the  joint  focal 
point,  where  the  amplitude  of  the  field  oscillates  at  the  difference  frequency  Aw.  The  region  of  interest  is 
placed  at  the  joint  focal  point  and  is  probed  point-by-point  by  raster  scanning.  The  sound  field  resulting 
from  object  vibrations  at  each  position  is  received  by  a  hydrophone  and  recorded.  The  recorded  signal 
at  one  or  more  difference  frequencies  is  used  to  form  an  image  of  the  object.  The  experiments  were 
conducted  in  a  water  tank.  The  transducer  center  frequency  was  3  MHz;  its  outer  diameter  was  45  mm; 
and  it  was  focused  at  70  mm.  The  difference  frequencies  used  in  each  experiment  are  mentioned  in  the 
corresponding  legends. 
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be  mapped  into  a  pictorial  format.  The 
spatial  resolution  of  the  resulting  image  is 
determined  by  the  region  in  which  signifi¬ 
cant  interference  between  the  ultrasound 
beams  occurs  and  is  of  the  order  of  a  few 
wavelengths  at  the  ultrasound  frequency. 

Experiments  were  conducted  in  a  water 
tank,  which  provided  good  ultrasonic  and 
acoustic  coupling  to  the  object  and  freedom 
of  movement  for  the  prototype  scanner 
mechanism  (Fig.  1).  The  two-element  con- 
focal  ultrasound  transducer  array  was  posi¬ 
tioned  such  that  the  beams  interfered  at  the 
selected  region  of  the  object.  Sound  pro¬ 
duced  by  vibrations  of  the  object  is  approx¬ 
imately  omnidirectional  because  of  the  small 
size  of  the  vibrating  portion  of  the  object 
compared  with  the  wavelength.  This  sound 
was  detected  by  a  submerged  hydrophone 
placed  near  the  object  within  the  water  tank. 

To  test  the  hypothesis  that  ultrasound- 
stimulated  acoustic  emission  is  sensitive  to 
object  mechanical  properties  and  to  show 
how  such  properties  can  be  quantitatively 
evaluated  by  this  method,  we  produced  an 
ultrasound-stimulated  vibro-acoustic  spec¬ 
trogram  of  a  tuning  fork  immersed  in  isopro¬ 
pyl  alcohol  at  two  different  temperatures. 
We  aimed  the  focal  point  of  the  confocal 
transducer  at  a  fixed  position  on  one  of  the 
tines.  The  shear  viscosity  of  alcohol  changes 
with  temperature,  causing  a  slight,  but  de¬ 
tectable,  shift  in  the  spectrogram  (Fig.  2). 
The  shear  viscosity  r\  of  a  liquid  is  deter¬ 
mined  by  measuring  the  resonant  frequency 
/R  and  the  bandwidth  S/R  of  a  tuning  fork 
immersed  in  this  liquid  (12) 


_  k/r/8/r  b/R0\2 

11  P  \  /r  /ro  / 


(3) 


where  /R0  and  8/R0  are  the  resonant  fre- 


Fig.  2.  Vibro-acoustic  spectrograms  of  a  tuning 
fork  immersed  in  isopropyl  alcohol  at  two  different 
temperatures.  A  point  on  a  tine  of  the  fork  was 
vibrated  with  the  use  of  the  system  shown  in  Fig. 
1 .  The  difference  frequency  was  swept  from  1 600 
to  2000  Hz.  A  change  of  temperature  from  1 5°  to 
30°C  decreases  the  shear  viscosity  of  the  alcohol, 
which,  in  turn,  changes  the  resonance  frequency 
and  bandwidth  of  the  tuning  fork. 


quency  and  bandwidth  measured  in  vacu¬ 
um,  respectively,  and  p  is  the  liquid  density. 
The  constant  k  is  determined  experimen¬ 
tally.  The  measured  values  for  /R  and  S/R 
were  1769  and  5.76  Hz  at  15°C  and  1765.7 
and  4-53  Hz  at  30°C.  The  viscosity  of  iso¬ 
propyl  alcohol  is  reported  to  be  2.89  cP  (1 
centipoise  =  1  mPa*s)  at  15°C  (13).  From 
this  value,  the  constant  k  was  calculated. 
The  shear  viscosity  at  30°C  was  found  using 
Eq.  3  to  be  T]  =  1.77  cP,  which  is  the  same 
as  the  published  data  (13). 

We  tested  the  ability  of  ultrasonically 
stimulated  vibro-acoustic  spcctrography  to 
image  the  frequency  response  of  different 
objects  with  identical  dr  by  scanning  three 
tuning  forks  with  different  resonant  frequen¬ 
cies.  A  color  acoustic  spectrogram  was  ob¬ 
tained  by  sweeping  the  frequency  of  the  ra¬ 
diation  force,  Aw,  in  a  range  covering  the 
resonant  frequencies  of  all  forks  at  each 
beam  position.  The  acoustic  emission  signal 
was  filtered  by  three  bandpass  filters  cen¬ 
tered  at  different  frequencies.  The  outputs  of 
these  filters  were  used  to  form  a  three-color 
composite  image  (Fig.  3).  The  forks  appear 


with  three  distinct  colors  because  each  fork 
responds  primarily  at  its  own  resonant  fre¬ 
quency.  Because  the  forks  were  made  from 
identical  materials,  other  ultrasound  imaging 
methods  would  not  be  capable  of  distin¬ 
guishing  these  objects. 

We  tested  the  hypothesis  that  the  ampli¬ 
tude  of  acoustic  emission  at  a  single  frequen¬ 
cy  can  be  used  to  detect  small,  highly  reflec¬ 
tive  isolated  objects.  We  scanned  a  380-fim- 
diameter  glass  bead  placed  on  a  thin  latex 
sheet  and  recorded  the  amplitude  of  acoustic 
emission  (Fig.  4).  The  latex  sheet  produces 
only  a  small  change  in  the  incident  energy 
because  it  is  almost  transparent  to  the  ultra¬ 
sound  beam.  This  experiment  demonstrated 
the  ability  of  the  method  to  detect  isolated 
regions  of  hardness  with  respect  to  a  soft 
background. 

To  test  the  feasibility  of  using  the  tech¬ 
nique  to  image  mechanical  properties  of 
tissues,  we  measured  the  phase  and  ampli¬ 
tude  of  acoustic  emission  from  calcified  and 
noncalcified  excised  human  iliac  arteries. 
The  arteries  were  scanned  in  a  plane  per¬ 
pendicular  to  the  ultrasound  beam  axis. 
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Fig.  3.  Vibro-acoustic 
spectrogram  of  tuning 
forks.  (A)  Three  tuning 
forks  made  from  identical 
material  with  identical  tine 
cross  sections  (lengths 
are  different).  Resonance 
frequencies  in  water  are 
407  Hz  (right),  809  Hz 
(middle),  and  1709  Hz 
(left).  The  forks  were 
scanned  in  a  water  tank 
with  the  use  of  the  sys¬ 
tem  shown  in  Fig.1 .  The 
scanning  plane  covers 

the  front  tines  at  the  bottom  part  of  the  forks.  At  each  position,  the  difference  frequency  was  swept  from 
250  to  2250  Hz.  The  ultrasound-stimulated  acoustic  emission  was  detected  with  the  hydrophone  and 
filtered  by  three  overlapping  bandpass  filters  with  frequencies  centered  at  500,  1000,  and  1500  kHz, 
respectively.  (B)  Color  acoustic  spectrogram  of  the  forks.  The  outputs  of  the  bandpass  filters  were  used  to 
produce  the  red,  green,  and  blue  image  components.  This  image  displays  two  characteristics  of  the 
object:  shape  and  frequency  response.  The  color  associated  with  each  fork  indicates  its  resonance 
frequency,  which  can  be  deduced  from  the  frequency  scale. 
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Fig.  4.  Image  of  ultrasonically  stimulated  acoustic 
emission  amplitude  from  a  380-fjim-diameter 
glass  bead  placed  on  a  thin  latex  sheet.  The  latex 
surface  was  scanned  at  0.2-mm  increments  in 
each  direction .  The  difference  frequency  was  fixed 
at  7.3  kHz.  The  amplitude  of  the  acoustic  emission 
of  the  bead  (in  relative  units)  is  shown  in  gray  < 

scale.  The  latex  sheet  is  almost  transparent  to  the  * 

imaging  system.  The  glass  bead,  however,  pre¬ 
sents  a  large  acoustic  impedance  discontinuity, 
resulting  in  significant  oscillatory  radiation  force. 

The  magnitude  of  the  radiation  force  gives  rise  to 
large-amplitude  acoustic  emission,  thus  yielding  a 
high-contrast  image  of  the  bead.  The  image  re¬ 
veals  a  system  resolution  of  about  700  p,m.  The 
confocal  beam  geometry  leads  to  a  negative-amplitude  ring  around  the  bright,  positive  central  spot.  This 
effect  produces  edge  enhancement,  as  seen  on  artery  walls  in  Fig.  5. 
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Calcifications  within  the  arteries  produced 
distinctive  amplitude  and  phase  values 
when  compared  to  the  normal  arterial  walls 
(Fig.  5).  The  phase  of  the  oscillation  of 
driven  mechanical  systems  relative  to  the 
driving  force  depends  on  the  ratio  of  mass 
to  stiffness  (14).  Calcified  regions  of  the 
diseased  artery,  identified  by  an  x-ray  of  the 
sample,  produced  phase  shifts  in  acoustic 
emission  completely  different  from  that  of 
the  noncalcified,  and  thus  softer  or  less 
dense,  regions.  The  amplitude  images  are 
highly  detailed  and  exhibit  variations  in 
acoustic  emission  from  both  calcified  and 
uncalcified  regions  of  the  diseased  artery. 
These  differences  are  caused  by  variations 
in  the  product  of  the  reflection  properties  dT 
and  the  effective  mechanical  vibration  ad¬ 
mittance  properties  Q(Aco)  of  the  tissue. 
Thus,  vibro-acoustic  spectrography  is  simi¬ 
lar  to  conventional  pulse-echo  ultrasound 
imaging,  which  is  sensitive  to  the  ultrasonic 
parameters  of  the  object,  but  has  the  advan¬ 
tage  of  also  being  sensitive  to  the  mechan¬ 
ical  admittance  Q(Aco)  at  low  frequencies. 

Motion  induced  by  ultrasound  and  mea¬ 
sured  with  ultrasound  pulse  echo  has  been 
used  previously  to  study  “hardness”  (15). 
However,  the  sensitivity  of  ultrasound  pulse 
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Fig.  5.  Vibro-acoustic  spectrography  of  excised 
human  iliac  arteries.  (A)  X-ray  image  of  normal  (left) 
and  calcified  (right)  excised  human  iliac  arteries 
obtained  from  a  35-year-old  woman  and  a  67- 
year-old  man,  respectively.  Bright  areas  indicate 
calcifications.  (B)  Vibro-acoustic  spectrogram 
amplitude  image  at  a  fixed  difference  frequency  of 
6  kHz.  Calcification  details  appear  bright,  whereas 
the  arterial  walls  are  dim.  (C)  Phase  image.  Calci¬ 
fied  regions  produce  acoustic  emission  of  differ¬ 
ent  phase  with  respect  to  regions  of  the  tissue 
having  little  calcification,  as  indicated  by  the  x-ray. 


echo  to  motion  at  common  ultrasound  fre¬ 
quencies  is  limited  to  several  micrometers. 
The  advantage  of  ultrasound-stimulated  vi¬ 
bro-acoustic  emission  is  its  high  displace¬ 
ment  sensitivity.  Cyclic  displacement  of 
100  nm  at  10  kHz  produces  an  acoustic 
intensity  of  about  3.0  X  10“ 3  W/cm2.  Hy¬ 
drophones  similar  to  the  one  used  in  these 
experiments  are  sensitive  to  as  little  as 
10“ 15  W/cm2  and,  therefore,  can  detect 
very  small  cyclic  displacements.  For  exam¬ 
ple,  the  hydrophone  detected  an  acoustic 
pressure  of  about  15  X  10“  3  Pa  at  a  distance 
of  5  cm  from  the  glass  bead  shown  in  Fig.  4- 
Under  the  assumption  of  isotropic  vibra¬ 
tion,  this  pressure  would  be  produced  by  a 
similarly  sized  sphere  vibrating  with  a  dis¬ 
placement  amplitude  of  about  6  nm.  The 
method  will  be  more  sensitive  for  higher 
frequency  sound  because  acoustic  power  is 
proportional  to  the  square  of  frequency  for 
constant  displacement  amplitude.  The 
practical  upper  limits  for  the  difference  fre¬ 
quency  produced  with  modern  ultrasound 
transducers  is  about  equal  to  their  band¬ 
width.  For  modern  transducers,  this  value  is 
80%  or  more  of  the  central  frequency  of  the 
transducer.  For  experiments  like  those  we 
conducted,  emission  frequency  well  in  ex¬ 
cess  of  1  MHz  could  be  produced.  The  lower 
limit  on  the  frequency  of  radiation  pressure 
is  zero,  that  is,  static  pressure. 

Ultrasound-stimulated  vibro-acoustic 
spectrography  has  potential  applications  in 
at  least  two  general  areas.  The  first  is  non¬ 
destructive  evaluation  of  materials,  where 
material  characteristics  and  structural  flaws 
can  be  identified  by  measuring  changes  in 
the  mechanical  response  to  vibration  at  a 
point.  The  object  under  test  could  be  re¬ 
motely  vibrated,  for  instance,  by  beams 
propagating  and  interfering  in  either  water 
or  air,  or  beams  propagating  within  the 
object  could  be  used  to  produce  acoustic 
emission  from  flaws.  For  medical  imaging 
and  detection,  the  technique  appears  par¬ 
ticularly  suitable  for  noninvasive  detection 
of  hard  tissue  inclusions,  such  as  the  imag¬ 
ing  of  arteries  with  calcification,  detection 
of  breast  microcalcifications,  visualization 
of  hard  tumors,  and  detection  of  foreign 
objects.  The  stiffness  of  soft  tissues  is  re¬ 
lated  to  their  composition  (for  example, 
relative  values  of  fibrotic  content),  and  its 
change  is  often  related  to  pathology  or 
therapy.  In  conventional  palpation,  phy¬ 
sicians  estimate  tissue  stiffness  by  feeling 
with  the  fingers.  Because  changes  of  stiff¬ 
ness  alter  the  vibration  frequency  response 
or  damping  of  tissue,  the  present  method 
can  potentially  provide  a  noninvasive,  re¬ 
mote,  high-resolution  “palpation”  tech¬ 
nique  that  can  reach  small  abnormalities 
that  are  otherwise  untouchable  by  con¬ 
ventional  methods. 
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The  long-term  time  average  is  obtained  when  T  -»<*. 
To  compute  the  short-term  time  average  of  the 
acoustic  energy  density  relevant  to  acoustic  emis¬ 
sion  at  Au>,  we  choose  T  longer  than  the  ultrasound 
wave  period  but  much  shorter  than  the  acoustic 
wave  period,  that  is,  2'tt/w2  <sc T  <sc2tt/A(d. 
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Abstract:  We  describe  theoretical  principles  of  a  new  imaging  modality  that  utilizes  the 
acoustic  response  of  an  object  to  a  highly  localized  dynamic  radiation  force  of  an  ultrasound 
field.  In  this  method,  named  Ultrasound-Stimulated  Vibro-Acoustography  (USVA),  ultrasound 
is  employed  to  exert  a  low-frequency  (in  kHz  range)  force  on  the  object.  In  response,  a  portion 
of  the  object  vibrates  sinusoidally  in  a  pattern  determined  by  its  viscoelastic  properties.  The 
acoustic  emission  field  resulting  from  object  vibration  is  detected  and  used  to  form  an  image 
that  represents  both  the  ultrasonic  and  low-frequency  (kHz  range)  mechanical  characteristics  of 
the  object.  We  report  the  relation  between  the  emitted  acoustic  field  and  the  incident  ultrasonic 
pressure  field  in  terms  of  object  parameters.  Also,  we  present  the  point  spread  function  of  the 
imaging  system.  The  experimental  images  in  this  report  have  a  resolution  of  about  700  //m,  high 
contrast,  and  high  signal-to-noise  ratio.  USVA  is  sensitive  enough  to  detect  object  motions  on 
the  order  of  nanometers.  Possible  applications  include  medical  imaging  and  material  evaluation. 

1  INTRODUCTION 

The  study  of  objects  in  terms  of  their  mechanical  response  to  external  forces  is  of  considerable 
interest  in  material  science  and  medical  diagnosis.  Elastic  constants  are  closely  connected  to  the 
thermodynamic  properties  of  materials,  and  can  be  related  to  a  wide  range  of  physical  parameters. 
Elastic  constants,  like  spring  constants,  can  be  determined  by  measuring  deformation  in  response 
to  an  applied  force.  Although  a  static  force  can  be  used  for  this  purpose,  using  a  dynamic  force 
is  preferred  if  one  is  interested  in  measuring  the  dynamic  characteristics  of  the  material  [1]. 

Changes  in  elasticity  of  soft  tissues  are  often  related  to  pathology.  Palpation  is  a  traditional 
example  of  estimating  mechanical  parameters  for  tissue  characterization,  where  a  static  force  is 
applied  and  a  crude  estimation  of  the  tissue  elasticity  is  obtained  through  the  sense  of  touch. 
In  palpation,  force  is  exerted  on  the  body  surface  and  the  result  is  an  accumulative  response 
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of  all  the  tissues  below.  Physicians  can  sense  abnormalities  if  the  response  to  palpation  of  the 
suspicious  tissue  is  sufficiently  different  from  that  of  normal  tissue.  However,  if  the  abnormality 
lies  deep  in  the  body,  or  if  it  is  too  small  to  be  resolved  by  touch,  then  the  palpation  method  fails. 

Elasticity  imaging,  a  subject  extensively  investigated  in  recent  years,  is  a  quantitative  method 
that  measures  the  mechanical  properties  of  tissue.  The  general  approach  is  to  measure  tissue 
motion  caused  by  an  external  (or,  in  some  methods,  internal)  force/displacement  and  use  it  to 
reconstruct  the  elastic  parameters  of  the  tissue.  Some  investigators  have  used  static  force  to 
compress  the  tissue  and  measured  the  resulting  strain  by  ultrasound  [2,  3].  Others  have  used 
external  mechanical  vibrators  to  vibrate  the  tissue  and  detected  the  resulting  displacement  in 
tissue  by  Doppler  ultrasound  [4-7].  For  a  review  of  elasticity  imaging  methods  refer  to  [8].  A 
recently  developed  method  employs  an  actuator  to  vibrate  the  body  surface  and  then  measures 
the  strain  waves  with  phase  sensitive  Magnetic  Resonance  Imaging  [9]. 

Most  of  the  elasticity  imaging  methods  are  based  on  an  external  source  of  force  resulting 
in  a  spatially  wide  stress  field  distribution.  This  requires  the  stress  field  to  pass  through  the 
superficial  portion  of  an  object  before  reaching  the  interior  part.  Analysis  of  the  object  response 
can  be  complicated  because  the  stress  field  pattern  changes,  often  unpredictably,  at  different 
depths  before  it  reaches  the  region  of  interest  within  the  object.  An  alternative  strategy  is  to 
apply  a  localized  stress  directly  in  the  region  of  interest.  One  way  to  accomplish  this  is  to  use 
the  radiation  pressure  of  an  ultrasound  source(s).  Based  on  this  strategy,  Sugimoto,  et  al.  [10] 
presented  a  method  to  measure  tissue  hardness  using  the  radiation  force  of  a  single  focused 
ultrasound  beam.  In  this  method,  impulsive  radiation  force  was  used  to  generate  localized 
deformation  of  the  tissue.  Resulting  transient  deformation  was  measured  as  a  function  of  time 
by  an  ultrasound  Doppler  method.  Radiation  force  has  also  been  used  to  generate  shear  elastic 
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waves  in  tissues  [11]. 

In  this  paper,  we  describe  the  principles  of  an  imaging  technique  that  produces  a  map  of  the 
mechanical  response  of  an  object  to  a  force  applied  at  each  point.  The  method  utilizes  ultrasound 
radiation  force  to  remotely  exert  a  localized  oscillating  stress  field  at  a  desired  frequency  within 
(or  on  the  surface  of)  an  object.  In  response  to  this  force,  a  part  of  the  object  vibrates.  The  size 
of  this  part  and  the  motion  pattern  depends  on  object  viscoelastic  characteristics.  The  acoustic 
field  resulting  from  object  vibration,  which  we  refer  to  as  “acoustic  emissions,”1  is  detected  by 
a  sensitive  hydrophone  and  used  to  form  the  image  of  the  object.  This  method  benefits  from 
the  high  spatial  definition  of  ultrasound  radiation  force  and  high  motion  detection  sensitivity 
offered  by  the  hydrophone.  We  call  this  technique  “ultrasound-stimulated  vibro-acoustography” 
(USVA).  Some  general  aspects  of  this  method  including  some  experimental  results  have  been 
outlined  by  the  authors  in  [12].  Here,  we  present  the  theoretical  foundations  of  USVA. 

2  METHODS 

Our  aim  is  to  image  an  object  based  on  its  mechanical  characteristics.  This  is  achieved  by 
vibrating  the  object  by  applying  a  highly  localized  oscillating  force  to  each  point  of  the  object. 
The  localized  force  is  produced  by  modulating  the  intensity,  and  thereby  the  radiation  force,  of 
the  ultrasound  at  low  frequencies  (normally  in  kHz  range).  The  resulting  sound  emitted  by  the 
object  is  a  function  of  object  mechanical  characteristics  and  the  location  of  the  excitation  point. 
The  image  is  produced  by  mapping  the  amplitude  or  phase  of  this  sound,  which  is  detected  by  a 
sensitive  hydrophone,  versus  position.  Figure  1  illustrates  this  method.  In  the  following  section 
we  describe  the  relationship  between  the  USVA  image  and  the  properties  of  the  object. 

1  The  term  “acoustic  emission”  is  used  to  describe  the  acoustic  field  in  response  to  a  cyclic  vibration  of  the  object.  Similar  terminology  is 
also  used  in  the  field  of  nondestructive  testing  of  materials  and  in  opto-acoustic  imaging  to  describe  a  different  phenomenon,  usually  the 
acoustic  field  resulting  from  structural  deformation. 
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3  THEORY 


Generation  of  a  dynamic  radiation  force  on  a  target 

The  acoustic  radiation  force  is  the  time  average  force  exerted  by  an  acoustic  field  on  an 
object.  This  force  is  an  example  of  a  universal  phenomenon  in  any  wave  motion  that  introduces 
some  type  of  unidirectional  force  on  absorbing  or  reflecting  targets  in  the  wave  path.  Radiation 
force  is  produced  by  a  change  in  the  energy  density  of  an  incident  acoustic  field.  For  a  review 
of  this  phenomenon  refer  to  [13].  Consider  a  plane  ultrasound  beam  interacting  with  a  planar 
object  of  zero  thickness  and  arbitrary  shape  and  boundary  impedance  that  scatters  and  absorbs. 
The  radiation  force  vector,  F,  arising  from  this  interaction  has  a  component  in  the  beam  direction 
and  another  transverse  to  it.  The  magnitude  of  this  force,  is  proportional  to  the  average  energy 
density  of  the  incident  wave  ( E )  at  the  object,  where  ()  represents  the  time  average,  and  S,  the 
area  of  the  projected  portion  of  the  object  [14] 


F  =  d  rS(E), 


(1) 


where  dr  is  the  is  the  vector  drag  coefficient  with  a  component  in  the  incident  beam  direction 
and  another  transverse  to  it.  The  coefficient  dr  is  defined  per  unit  incident  energy  density  and 
unit  projected  area.  For  a  planar  object,  dr  is  numerically  equal  to  the  force  on  the  object. 
Physically,  the  drag  coefficient  represents  the  scattering  and  absorbing  properties  of  the  object 
and  is  given  by  [14] 


(ria  +  IL, 


/ 


7  cos  asdS 


7  sin  asdS, 


(2) 


where  p  and  q  are  the  unit  vectors  in  the  beam  direction  and  normal  to  it,  respectively.  The 
quantities  ITa  and  IIg  are  the  total  absorbed  and  scattered  powers,  respectively,  and  7  is  the 
scattered  intensity,  all  expressed  per  unit  incident  intensity.  Also,  as  is  the  angle  between  the 
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incident  and  the  scattered  intensity,  and  dS  is  the  area  element.  The  drag  coefficient  can  also  be 
interpreted  as  the  ratio  of  the  radiation  force  magnitude  on  a  given  object  to  the  corresponding 
value  if  the  object  were  replaced  by  a  totally  absorbing  object  of  similar  size.  This  is  because 
|  dr  |=1  for  a  totally  absorbing  object.  This  coefficient  can  be  determined  for  objects  of 
different  shapes  and  sizes.  For  simplicity,  we  assume  a  planar  objected  normal  to  the  beam  axis. 
In  this  case,  the  transverse  component  vanishes,  thus,  the  drag  coefficient  (force)  will  have  only 
a  component  normal  to  the  target  surface  which  we  denote  by  scalar  dr  (F).  Values  of  dr  for 
spheres,  in  terms  of  the  diameter  and  the  wavelength,  are  given  in  [14]. 

To  produce  a  dynamic  radiation  force,  one  can  use  an  amplitude  modulated  beam  [15]. 
Consider  an  amplitude  modulated  incident  (ultrasonic)  pressure  field,  p(t),  as 

p(t)  =  PWo  cos  (Awt/2)  cos  Dot. ,  (3) 

where  PWo,  Aw/2,  and  w o  are  the  pressure  amplitude,  modulating  frequency,  and  the  center 

frequency,  respectively.  In  our  analysis  and  experiments,  we  assume  that  the  condition  Auj  « 

wo  holds.  In  such  a  case,  the  energy  density  of  the  incident  field  has  slow  variations  in 

time.  To  discriminate  the  slow  time  variations  of  a  function,  let  us  define  the  short-term 

time  average  of  an  arbitrary  function  £(t)  over  the  interval  of  T  seconds  at  time  instance 

t,  as  (£( t))T  =  y  f  £( r)dr ,  which  is  a  function  of  t.  The  long-term  time  average  (or 
t—T/2 

simply,  the  time  average)  is  obtained  by  setting  T  -+  oo.  To  compute  the  short-term  time 
average  of  the  acoustic  energy  density  relevant  to  field  variations  at  Aw/2,  we  choose  T 
longer  than  the  ultrasound  wave  period  but  much  shorter  than  the  modulation  period,  that  is 
2n/u>o  «  T  «  A-k/Aw.  Under  this  condition,  the  short-term  time  average  of  p2(t)  is 
\P  (t))?  =  +  cosA wt).  The  energy  density  is  given  by  p2(t)/pc 2  where  p  and  c  are 

the  density  and  propagation  speed  in  the  medium  [16].  We  are  interested  in  the  time-varying 
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component  of  the  short-term  time  average  of  the  energy  density.  Denoting  this  component  by 
eAui(t),  we  can  write:  eAuJ(t)  =  cos  Amt.  This  component  of  the  energy  density  produces 
a  time  varying  radiation  force  on  the  target  (Eq.  1)  at  frequency  Au>.  The  amplitude  of  this 
force,  FAuj,  is 

Fauj  =  4pc*Sdr' 

This  equation  states  that  the  time-varying  force  amplitude  is  proportional  to  the  square  of  incident 
ultrasound  pressure,  or  equivalently,  to  the  incident  power.  If  the  object  moves  in  response  to 
this  force,  then  the  high-frequency  ultrasound  energy  would  convert  to  low-frequency  mechanical 
energy. 

Acoustic  emission  from  a  target  due  to  a  dynamic  force 

The  radiation  force  FAlj  vibrates  the  target  object  at  frequency  Au /.  Object  vibration  results 
in  an  acoustic  field  in  the  medium  (acoustic  emission).  This  field  is  related  to  object  shape, 
size,  and  viscoelastic  properties.  To  present  a  conclusive  analysis  of  this  relationship,  we  have 
to  assume  an  object  with  specific  characteristics.  Here,  we  assume  that  the  vibrating  object  has 
a  circular  cross  section  of  radius  b  and  uniformly  vibrates  back  and  forth  like  a  piston.  This 
choice  allows  us  to  illustrate  the  concept  in  a  simple  form.  We  also  consider  an  area  S  <  7 rb2 
of  the  piston  surface  to  be  projected  normally  by  the  beam.  Similar  solutions  can  be  carried 
out  for  other  objects.2 

The  steady  state  normal  velocity  amplitude  of  a  piston,  UAu},  due  to  a  harmonic  force  FAu) 
at  frequency  Ao>,  can  be  described  in  terms  of  the  mechanical  impedance  ZAu), 

TT  F 

_  Uau  ~  zZ’  (5) 

2  The  theory  can  be  extended  to  include  arbitrary  vibrating-part  shapes  and  nonuniform  displacement  of  the  object.  Nonuniform  displacement 

would  be  an  important  issue  when  the  vibration  wavelength  in  the  object  material  is  smaller  than  2b. 
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where  Z aw  —  Z'm  +  ZT  is  comprised  of  the  mechanical  impedance  of  the  object  in  vacuum 
Z'm,  and  the  radiation  impedance  of  the  object  Zr,  all  defined  at  Aw.  Modeling  the  object  as  a 
mass-spring  system,  Z'm  can  be  written  in  terms  of  Auj  as  [16,  17] 

z'm  =  R'm-  3  ,  (6) 

where  m,  R'm,  and  K'  are  the  mass,  mechanical  resistance,  and  the  spring  constants  of  the  object, 
respectively.  The  radiation  impedance  of  the  piston  can  be  written  as  [17] 

Zr  =  *b\Rr-jXr),  (7) 


where 


Rr  =  pc\l  -  ——rJ\( — - — )1, 
H  L  A ub  1\2Au;bJy 


Xr  = 


'/sin( 


2A  ujb 


cos  a  )  sin2  ada, 


where  Ji(-)  is  the  first-order  Bessel  function  of  the  first  kind.  In  many  applications  of  our 
interest,  the  wavelength  is  much  greater  than  the  object  size,  hence,  — b  — ►  0.  In  such  cases 
Zr  assumes  a  simpler  form: 

7  ,3  .  / A u>b  .  8  \ 

Z,  =  *bPA^—-3-).  (10) 

The  mechanical  impedance  of  the  piston  object  can  now  be  written  as 
Zalj  =  ( R'm  +  V b2Rr)  -  j  Au  -  ^ -  -I-  7 rb2Xr^j 

Once  we  calculate  U^,  we  can  calculate  the  pressure  field  it  produces  in  the  medium.  We 
assume  that  the  acoustic  emission  signal  propagates  in  a  free  and  homogenous  medium.  The 
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farfield  acoustic  pressure  due  to  a  piston  source  of  radius  b  set  in  a  planar  boundary  of  infinite 
extent  is  given  by  [17], 


Paw  =  -jAup 


exp  (jAtvl/c) 
47 rZ 


c  sin^)  cost? 


A  ujb 
c 


sin 'd 


x 


COS  7?  +  Pb 


(2Trb2UAuJ),  (12) 


where  l  is  the  distance  from  the  observation  point  to  the  center  of  the  piston,  i?  is  the  angle 
between  this  line  and  the  piston  axis,  and  Pb  is  the  specific  acoustic  admittance  of  the  boundary 
surface.3  The  factor  of  two  comes  from  the  presence  of  the  boundary  wall.  It  would  be  replaced 


by  unity  if  the  boundary  wall  were  not  present  [16].  The  acoustic  emission  field  resulting  from 
object  vibration  can  be  written  in  terms  of  the  incident  ultrasound  pressure  by  combining  Eqs.  4, 
5,  and  12,  as 


x 


exp  (  jAtul/c) 
47 rl 

1 


27i(^^sini?)  cos  t? 


sin 


(R'm  +  7T b2Rr)  -  jlmAuj  -  +  tt b2Xr) 


cost?  +  Pb 

{2*b2)r*,sdr, 


(13) 


For  wavelengths  long  compared  to  the  object  size,  i.e.,  when  bAu/c  ->  0,  the  term  in  the 
first  brace  approaches  a  constant,  hence  we  may  consider  the  contents  of  the  first  brace  to  be 
an  object  independent  function  (the  specific  acoustic  admittance  Pb  relates  to  the  surrounding 
boundary  surface).  Under  these  conditions,  the  first  brace  in  the  above  equation  represents  the 
effect  of  the  medium  on  the  acoustic  emission  field,  which  we  may  call  the  medium  transfer 
function,  and  denote  it  by 


Hau,(1) 


Au  exp  (jAojl/c) 
c2  47 rl 


2ji(A^sin  i9)  cos 


A  ujb 
c 


sin$  cost 9  +  Pb 


(14) 


The  specific  acoustic  admittance  is  0 b  —  *  where  Zg,  the  acoustic  impedance  of  the  boundary,  represents  the  ratio  between  the 

pressure  and  normal  fluid  velocity  at  a  point  on  the  surface. 
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The  second  bracket  in  Eq.  13  is  z~t  or  the  mechanical  admittance  of  the  object  at  the  frequency 
of  the  acoustic  emission  (Aw),  and  we  denote  it  by  YAu).  It  is  convenient  to  combine  this  term 
with  the  next  term  (2ttZ>2)  in  Eq.  13,  as  QAuJ  =  2tt  b2YAuJ  =  2tt b2/ZAuJ,  which  is  the  total  acoustic 
outflow4  by  the  object  per  unit  force.  Function  QA(jJ  represents  the  object  characteristics  at  the 
acoustic  frequency.  We  may  thus  rewrite  Eq.  13  in  a  more  compact  form  as 

PAw  =  HAuJ{l)QAwP2QSdT.  (15) 

Equation  15  indicates  that  the  acoustic  emission  pressure  is  proportional  to:  1)  the  square  of 
ultrasound  pressure  Pwo,  2)  the  ultrasound  characteristics  of  the  object,  dr,  in  the  projected  area 
S,  3)  the  acoustic  outflow  by  this  object,  QAjJ,  representing  the  object  size  b  and  its  mechanical 
admittance  at  the  acoustic  frequency,  YAu},  and  4)  the  transfer  function  of  the  medium  at  the 
acoustic  frequency,  HAu)(l).  The  above  equation  illustrates  the  basic  nonlinear  relationship 
between  the  ultrasound  and  acoustic  emission  pressure  amplitudes.  Note  that  neither  the  medium 
nor  the  object  need  to  be  nonlinear  for  this  relationship  to  hold.  It  is  interesting  to  note  that  the 
projection  area  S,  and  the  vibrating  area  ixb2  play  different  roles.  The  projection  area  determines 
the  extent  of  the  force  applied  to  the  object  (Eq.  4).  The  vibrating  area,  however,  influences 
the  total  acoustic  outflow  in  the  medium  caused  by  object  vibration.  The  mechanism  of  object 
vibration  is  somewhat  analogous  to  that  of  a  loudspeaker,  where  the  electromotive  force  is 
exerted  to  a  small  area  of  the  membrane  (usually  at  the  center)  causing  the  entire  membrane 
surface  to  vibrate.  In  our  method,  the  size  of  the  vibrating  area  depends  on  the  object  structure. 
For  a  free  suspended  point  object,  smaller  than  the  beam  cross  section,  the  vibrating  area  would 
be  the  same  as  the  projection  area.  For  a  large  stiff  plate,  however,  the  vibrating  area  could 

Acoustic  outflow  is  the  volume  of  the  medium  (e.g.,  the  fluid)  in  front  of  the  object  surface  that  is  displaced  per  unit  time  due  to  object 

vibration. 
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be  much  larger  than  the  projected  area  (similar  to  a  loudspeaker).  In  some  cases  it  is  more 
convenient  to  have  the  acoustic  emission  field  in  terms  of  the  applied  force  F Referring  to 
Eq.  4,  we  can  rewrite  Eq.  15  as 

PAu>  =  4pc2HAuJ(l)QAu,FAuJ-  (16) 

Again  in  analogy  to  a  loudspeaker,  i?Ao>>QAuMan(l  #Aw(0  represent  the  electromotive  force, 
dynamic  characteristics  of  the  membrane,  and  the  propagation  medium  transfer  function. 

Beam  forming 

To  probe  an  object  with  the  dynamic  radiation  force  at  high-spatial  resolution,  it  is  ideal  to 
confine  the  dynamic  stress  field  to  a  very  small  region  in  three-dimensional  space.  We  may  define 
the  resolution  cell  of  the  system  as  the  volume  within  which  the  amplitude  of  the  modulated  field 
is  high  enough  to  produce  a  stress  field  on  a  target.  The  purpose  of  beam  forming  is  to  produce  a 
resolution  cell  as  small  as  possible.  An  amplituded  modulated  single-focused  beam  can  provide 
a  resolution  cell  that  is  small  in  diameter  but  long  in  the  depth  direction.  A  superior  strategy  that 
can  achieve  small  resolution  cell  in  all  dimensions  is  to  use  two  unmodulated  focused  beams 
at  slightly  different  frequencies,  and  allow  them  to  cross  each  other  at  their  focal  regions.  This 
is  accomplished  by  projecting  two  coaxial,  confocal,  continuous  wave  (CW)  ultrasound  beams 
on  the  object.  An  amplitude  modulated  field  is  produced  only  at  the  interference  region  of  the 
two  unmodulated  beams  around  their  focal  areas,  resulting  in  a  small  resolution  cell.  For  this 
purpose,  elements  of  a  two-element  spherically  focused  annular  array  (consisting  of  a  central 
disc  with  radius  ai  and  an  outer  ring  with  the  inner  radius  of  a'2  and  outer  radius  of  a^)  are 
excited  by  separate  CW  signals  at  frequencies  wi  =  wo  -  Aw/2  and  W2  =  wo  +  Aw/2.  We 
assume  that  the  beams  are  propagating  in  a  lossless  medium,  in  the  +z  direction  of  a  Cartesian 
coordinate  system  ( x,y,z ),  with  the  joint  focal  point  at  2  =  0.  The  resultant  pressure  field  on 
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the  2  =  0  plane  may  be  written  as 


p(t)  =  Pi(r)cos(ioit  +  +  P2(r)  cos(u2t,  +  ^(r)), 

where  r  =  ^/x2  +  y2  is  the  radial  distance.  The  amplitude  functions  are  [16,  18] 

Pi(r)  =  pcUov 
and 


(17) 


P2(r)  =  pcUo2 
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where  A =  27r/w,;,  i=],  2,  is  the  ultrasound  wavelength,  Uqi  is  the  particle  velocity  amplitude 
at  the  i-th  transducer  element  surface,  and  jinc(X)  =  Ji(2irX)/irX.  The  phase  functions, 
i’i  (r )  =  —  for  i= 1,  2,  are  conveniently  set  to  be  zero  at  the  origin. 

Now,  we  define  a  unit  point  target  at  position  (xo,ijq)  on  the  focal  plane  with  a  drag 
coefficient  distribution  as 

dr(x,y)  =  6(x-x0,y-yo),  (20) 

such  that  dr(x,y)dxdy  is  unity  at  the  (xo,y0)  and  zero  elsewhere.  This  equation  is  merely  used 
as  a  mathematical  model  because  dr  is  physically  finite.  In  this  case,  the  projected  area  can  be 
considered  to  be  S  =  dxdy.  We  replace  dTS  in  Eq.  1  with  dr(x,y)dxdy,  and  follow  the  steps 
similar  to  those  outlined  in  Eqs.  3  and  4  for  the  pressure  field  expressed  by  Eq.  17,  then  the 
complex  amplitude  of  the  normal  component  of  the  force  on  the  unit  point  target  can  be  found  as 
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where  the  arguments  xq  and  yo  are  added  to  denote  the  position  of  the  point  target,  and 
ro  =  \Jxl  +  Vo-  Equation  21  describes  the  spatial  distribution  of  the  force  (the  stress  field).  This 
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equation  shows  that  the  stress  field  is  confined  to  the  regions  near  the  beam  axis  (r0  =  0),  and 
decays  as  the  radial  distance  r0  increases.  The  lateral  extent  of  the  stress  field,  and  hence  the 
resolution  cell  diameter,  would  be  smaller  at  higher  ultrasound  frequencies  (smaller  Ai  and  A2). 
One  can  calculate  the  total  force  on  an  arbitrary  object  by  integrating  the  force  over  the  projected 
area.  The  axial  extent  of  the  resolution  cell  (depth  resolution,  or  the  depth  of  field)  can  be 
determined  by  calculating  the  force  FAoJ  as  a  function  of  the  depth  variable,  in  a  fashion  as 
outlined  in  Eqs.  17  to  21.  For  conciseness,  we  will  only  present  the  measured  values  for  the 
depth  resolution  in  the  Results  section. 

Loss  in  the  propagation  path  would  attenuate  both  ultrasound  beams,  thus  less  radiation 
force  would  be  generated  by  the  remaining  ultrasound  energy.  In  the  case  of  soft  tissues,  the 
force  attenuation  factor  is  A(zq)  =  exp  [<x.zq{uj\  +  w2)],  where  a  is  the  attenuation  coefficient  of 
the  tissue.  Energy  loss  in  the  medium  would  also  result  in  generation  of  a  separate  radiation 
stress  on  the  medium  along  the  ultrasound  paths.  However,  since  the  two  beams  propagate  along 
separate  paths  in  the  CW  form,  they  exert  mainly  steady  radiation  stresses  to  the  medium,  which 
does  not  cause  object  or  medium  vibrations.  Dynamic  radiation  force  is  produced  only  in  the 
interference  region  around  the  focal  area.  This  is  another  advantage  of  using  two  unmodulated 
beams  over  a  modulated  single  beam. 

Image  formation 

To  produce  an  image  we  scan  the  object  in  a  plane  and  record  the  complex  amplitude  of 
the  acoustic  emission,  at  different  positions.  In  this  process,  we  keep  Aw  fixed.  For 
transverse  view  images,  the  scan  plane  is  the  focal  plane  (x-y).  Alternatively,  for  the  parallel 
view  the  scan  plane  is  the  x  —  z  plane.  In  the  conventional  ultrasound  imaging  context,  these 
two  views  are  called  the  C-scan  and  B-scan,  respectively.  Our  main  focus  here  is  the  transverse 
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view  imaging.  In  this  case,  the  acoustic  emission  data  obtained  by  vibrating  the  object  at  point 
(x,  y )  are  assigned  to  the  corresponding  point  (z,  y)  in  the  image. 

Before  defining  the  image,  we  need  to  define  the  function  that  represents  the  object.  Referring 
to  Eq.  15,  the  terms  that  are  object  dependent  are  the  drag  coefficient  dr  and  the  function  Qau 
(assuming  that  HAu)(l )  is  object  independent).  The  object  function  g(x,  y)  is  thus  defined  as 
the  spatial  distribution  of  these  terms, 


g(x,  y)  =  QAu,(z,  y)dr(x,  y).  (22) 

Variables  x  and  y  are  added  to  denote  the  dependency  of  d.T  and  on  position.  In  particular, 
QAu>(x,y)  implies  the  total  acoustic  outflow  by  the  object  when  the  force  is  applied  at  point 
(*,  y )• 

Commonly,  an  imaging  system  is  studied  through  its  point  spread  function  (PSF),  which  is 
defined  as  the  image  of  a  point  object.  To  determine  the  PSF  of  our  system,  we  consider  a  unit 
point  target  at  the  origin  with  unit  mechanical  response,  Q&w{x,y)  =  1.  Hence,  referring  to 
Eqs.  20  and  22,  we  can  write  g(x,y)  =  6(x,y).  To  obtain  the  PSF,  we  move  this  point  object 
to  every  possible  position  (zo,  yo)  on  the  z  =  0  plane  and  form  the  image  using  the  resulting 
acoustic  emission  field,  Pau)(^0,Vo)-  Since  zq  and  yo  are  now  being  treated  as  variables,  we 
may  replace  them  by  variables  z  and  y,  respectively.  We  define  the  normalized  PSF  of  the 
coherent  imaging  system  as  the  complex  function 


h(x,y)  =  PAu(x,y)/PAu(0,0). 


(23) 


Division  by  Pau,(0,  0)  cancels  the  constant  multipliers.  Referring  to  Eqs.  16  and  21,  we  can  write 
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This  equation  illustrates  that  the  system  PSF  is  a  circularly  symmetric  function  with  the  peak  at 
the  origin  and  decaying  amplitude  with  increasing  the  radial  distance  r.  Amplitude  decays  faster 
for  higher  ultrasound  frequency.  This  function  will  be  discussed  further  in  the  next  section. 

4  EXPERIMENTS 

The  experimental  setup  is  shown  in  Fig.  2.  The  confocal  transducer  is  constructed  using  a 
spherical  piezoelectric  cap.  The  two  elements  are  constructed  by  dividing  the  back  electrode 
of  the  piezoelectric  wafer  into  a  central  disc  and  the  outer  ring,  such  that  the  elements  have 
identical  beam  axes  and  focal  lengths.  Radii  of  the  elements  are  o.i=14.8  mm,  a2= 22.5  mm,  and 
0.2=16.8  mm,  and  the  focal  distance  is  70  mm.  Transducer  elements  were  driven  by  two  stable 
RF  synthesizers  (HP  33120A  and  Analogic  2045)  at  frequencies  of  /0  -  A//2  and  /o  +  A//2, 
where  /o=  3  MHz,  and  the  value  of  A /  =  Au>/2n  is  stated  separately  for  each  experiment.  The 
object  was  placed  at  the  focal  plane  of  the  ultrasound  beams  in  a  water  tank.  Sound  produced 
by  the  object  vibration  was  detected  by  an  audio  hydrophone  (ITC  model  680,  sensitivity  -154 
dB  re  lV//iPa)  placed  within  the  water  tank.  The  received  signal  was  filtered  and  amplified 
by  a  programmable  filter  (Stanford  Research  Systems,  SR650)  to  reject  noise,  then  digitized  by 
a  12  bits/sample  digitizer  (HP,  E1429A)  at  a  rate  sufficiently  higher  than  the  Nyquist  rate  for 
the  particular  A /  used.  Data  were  recorded  on  a  computer  disc  (Sun  Sparc  computer).  For 
coherent  imaging,  which  requires  the  phase  information,  the  reference  signal  (i.e.,  cos  A  cot)  was 
obtained  by  electronic  down  mixing  of  the  two  driving  signals  and  was  recorded  along  with 
the  hydrophone  signal.  The  relative  phase  of  the  acoustic  emission  data  was  then  calculated  at 
each  point  by  discrete  Hilbert  transform.  We  conducted  two  experiments.  The  first  experiment 
is  designed  to  verify  the  relationship  between  the  acoustic  emission  pressure  and  the  ultrasound 
pressure  (Eq.  15),  and  the  second  is  designed  to  experimentally  measure  the  point  spread  function 


41 


stated  in  Eq.  24.  These  experiments  are  presented  to  prove  the  principles  of  the  method.  Further 
experiments,  illustrating  applications  of  the  method  for  evaluation  of  object  mechanical  properties 
and  tissue  imaging,  are  presented  in  [12]. 

5  RESULTS 

Acoustic  emission  versus  ultrasound  pressure 

Equation  15  states  that  the  acoustic  emission  field  amplitude  is  linearly  proportional  to  the 
square  of  the  incident  ultrasound  pressure,  or  equivalently,  to  the  incident  power.  To  test  this 
hypothesis,  a  calibrated  1  mm  diameter  ultrasound  needle  hydrophone,  with  its  tip  facing  the 
ultrasound  beam,  was  used  to  measure  the  ultrasound  field  at  the  focal  point.  (This  is  in  addition 
to  the  audio  hydrophone.)  The  tip  of  the  hydrophone  also  served  as  an  object  to  generate 
the  acoustic  emission.  Here,  A /  was  set  at  40  kHz,  and  the  radial  distance  from  the  tip  of 
the  needle  hydrophone  to  the  audio  hydrophone  was  about  65  mm.  The  result  is  shown  in 
Fig.  3.  The  slope  of  the  acoustic  emission  intensity  versus  ultrasonic  intensity  indicates  that  the 
intensity  of  the  acoustic  emission  field  is  proportional  to  the  square  of  the  ultrasound  intensity,  or 
equivalently,  the  acoustic  emission  pressure  amplitude  is  linearly  proportional  to  the  ultrasound 
power,  as  predicted  by  Eq.  15.  In  another  experiment,  a  450  /zm  diameter  glass  bead  was  used 
as  a  point  object.  In  this  case,  A /  was  set  at  either  7  or  40  kHz.  The  radial  distance  from 
the  glass  bead  to  the  audio  hydrophone  was  about  50  mm.  Again,  the  data  indicate  a  quadratic 
relationship  between  the  acoustic  emission  and  ultrasonic  intensities.  These  glass  bead  data  also 
show  that  increasing  the  frequency  increases  the  acoustic  intensity.  This  can  be  better  understood 
by  investigating  the  theoretical  model  presented  in  Eq.  13.  If  we  assume  that  the  mechanical 
admittance  of  the  object  is  approximately  constant  at  these  frequencies,  then  the  object  behaves 
almost  as  a  point  source  and  the  acoustic  pressure  amplitude  is  proportional  to  Aw.  Hence, 
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the  acoustic  emission  intensities  I^Wx  and  7ao>2  at  frequencies  Awi  and  Au/2,  respectively,  are 
related  by  /Aa>2//Awi  =  (Aw2/Awi)2.  Now,  letting  Awi  and  Au/2  correspond  to  7  and  40  kHz, 
respectively,  we  can  write  the  intensity  ratio  as:  [IaujJIa^B  =  201oS  (Au;2/Ao>i)  =  15  dB. 
The  mean  value  of  the  intensity  ratio  calculated  from  the  glass  bead  data  at  7  and  40  kHz  is  16 
dB,  which  is  in  close  agreement  with  the  theoretical  result. 

PSF  Measurement 

To  demonstrate  the  image  formation  process  and  support  the  theoretical  derivation  of  the 
PSF  (Eq.  24),  we  evaluated  this  function  experimentally.  For  this  purpose,  we  used  a  380  //m 
diameter  glass  bead  as  a  model  for  a  point  and  placed  it  on  a  thin  latex  sheet.  The  latex  sheet 
produces  only  a  small  change  in  the  incident  energy,  and  hence,  does  not  produce  significant 
radiation  force  or  acoustic  emission.  The  entire  object  was  placed  in  a  water  tank  and  the 
latex  sheet  surface  was  scanned  in  a  raster  format  at  0.2  mm  increments  in  either  direction  at 
Af-7.3  kHz.  The  amplitude  and  phase  of  the  acoustic  emission  signal  were  calculated  at  each 
point  relative  to  the  reference  signal  data.  The  phase  was  normalized  to  the  phase  value  at 
the  center  of  the  bead.  The  resulting  in-phase,  quadrature,  phase,  and  magnitude  images  are 
shown  in  Figs.  4a  through  d.  Transverse  image  resolution,  defined  as  the  -6dB  width  of  the  bead 
image,  is  approximately  700  fj.m  for  the  in-phase  image  in  either  dimension  (refer  to  Fig.  5).  To 
compare  the  experimental  results  with  those  of  the  theory,  we  calculated  the  profile  of  the  PSF, 
for  the  transducer  parameters  used  in  this  experiment,  according  to  Eq.  24.  Figure  5  shows  the 
theoretical  PSF  profile  and  the  glass  bead  profile  obtained  from  the  experiment  (in-phase  image 
profile,  shifted  to  center  at  zero).  This  figure  shows  excellent  agreement  between  theory  and 
experiment  for  amplitudes  above  20%  of  the  peak.  The  experimental  data  show  some  background 
offset  about  12%  of  the  peak.  We  believe  that  this  background  offset  is  caused  at  least  by  the 
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following  sources:  1)  acoustic  emission  by  the  latex  sheet,  2)  the  background  acoustic  noise  in 
the  experimental  setup  caused  by  equipment  fans  and  some  structural  and  building  vibrations,  3) 
nonlinearity  of  water  which  can  produce  a  nonlinear  mix  of  the  two  beams  even  in  the  absence 
of  the  object,  and  4)  streaming  as  a  result  of  energy  absorption  by  water  [14]  which  in  turn 
can  vibrate  the  latex  sheet  blocking  the  stream.  To  evaluate  the  depth  resolution  (or  the  slice 
thickness),  we  placed  the  glass  bead  on  the  beam  axis  and  scanned  it  in  the  z-direction  about  the 
focal  point.  The  depth  resolution,  defined  as  the  distance  between  the  points  where  the  amplitude 
of  the  acoustic  emission  field  drops  to  -6  dB  of  its  peak,  was  9  mm. 

6  DISCUSSION 
System  properties 

Spectral  Characteristics — In  general,  USVA  images  represent  object  characteristics  at  two 
ends  of  the  spectrum:  the  drag  coefficient  at  the  ultrasound  frequency,  and  the  mechanical 
admittance  at  the  low  acoustic  frequency  A/.  The  ultrasound  frequency  is  usually  set  at  a  value 
suitable  to  form  the  beam,  while  A /  can  vary  in  a  wide  range  depending  on  the  application.  If 
the  two  beams  are  produced  by  similar  ultrasound  transducer  elements,  then  the  practical  upper 
limit  for  A /  is  about  equal  to  the  transducer  bandwidth.  The  lower  limit  of  A /  is  zero. 

Sensitivity — System  sensitivity  in  detecting  very  small  displacements  is  an  important  practical 
issue,  especially  when  the  allowable  ultrasound  power  is  limited  (for  example,  in  medical 
imaging).  Motion  measurement  with  ultrasound  pulse  echo  has  been  used  previously  to  study 
“stiffness”  of  tissues  [10].  However,  the  sensitivity  of  ultrasound  pulse  echo  to  motion,  at 
common  medical  ultrasound  frequencies,  is  limited  to  several  micrometers.  An  advantage  of 
USVA  is  its  high  displacement  sensitivity.  Cyclic  displacement  of  100  nanometers  at  10  kHz 
produces  an  acoustic  intensity  of  about  3.0x1 0"3  watts/cm2.  Hydrophones  similar  to  the  one 


used  in  these  experiments  are  sensitive  to  as  little  as  10~15  walls/cm2,  and  therefore,  are  capable 
of  detecting  very  small  cyclic  displacements.  For  instance,  the  hydrophone  detected  an  acoustic 
pressure  of  about  15xlO-3  Pascals  at  a  distance  of  5  cm  from  the  glass  bead  shown  in  Fig.  4. 
Under  assumptions  of  isotropic  vibration,  this  pressure  would  be  produced  by  a  similar  sized 
sphere  vibrating  with  a  displacement  amplitude  of  about  6  nm.  Sensitivity  increases  at  higher 
frequencies  because  the  acoustic  emission  pressure  is  proportional  to  frequency  for  constant 
mechanical  admittance  (Eq.  13). 

Comparisons  with  pulse-echo  systems — Some  contrasting  features  of  USVA  with  respect  to 
the  conventional  ultrasound  pulse-echo  imaging  (B-mode  and  C-mode)  are: 

1.  A  pulse-echo  image  represents  object  microstructure  by  displaying  its  ultrasonic  reflectivity 
distribution.  The  acoustic  emission  signal  in  a  USVA  system  is  proportional  to  the  drag 
coefficient  dr(x,y),  which  is  a  local  ultrasound  parameter,  and  to  the  function  Q&u>{x,y), 
which  represents  the  bulk  response  of  the  object  at  acoustic  frequency  Aw.  Hence,  a  USVA 
image,  in  general,  represents  both  the  microstructure  and  macrostructure  of  the  object. 

2.  The  echo  signal  in  a  pulse-echo  system  is  a  linear  function  of  the  incident  ultrasound  pressure 
amplitude  and  the  amplitude  reflection  coefficient  of  the  object.  In  USVA,  however,  the 
acoustic  emission  signal  is  proportional  to  the  ultrasound  power  and  the  power  reflection 
coefficient  of  the  object. 

3.  Pulse-echo  systems  are  not  directly  sensitive  to  medium  absorption.  Absorption  is  indicated 
as  relative  changes  in  the  amplitude  of  the  echoes  resulting  from  the  scatterers  within  the 
medium.  In  a  USVA  system,  the  acoustic  emission  can  be  produced  directly  as  a  result  of 
energy  absorption  by  the  medium,  even  if  the  medium  is  homogeneous  (refer  to  Eq.  2). 

4.  Pulse-echo  systems  are  generally  broadband.  The  USVA  method  presented  here  is  basically 


a  narrowband  technique. 

5.  Pulse  echo  systems  achieve  high  depth  resolution  by  transmitting  short  wideband  pulses.  A 
USVA  system  gains  its  depth  resolution  by  tailoring  beam  geometry  to  limit  the  depth  of  the 
region  where  the  two  beams  interfere.  A  USVA  system  does  not  require  a  wide  bandwidth 
signal  to  achieve  a  high  depth  resolution. 

6.  The  data  of  USVA  images  are  acquired  one  point  at  a  time,  which  resembles  the  data 
acquisition  in  C-mode  pulse-echo  systems.  B-mode  pulse-echo  systems,  however,  require 
much  less  acquisition  time  because  the  data  are  collected  one  line  at  a  time. 

Applications 

The  method  presented  here  promises  applications  in  two  general  areas:  medical  imaging 
and  material  evaluation. 

Medical  Applications — USVA  can  be  used  to  image  tissues  and  evaluate  their  mechanical 
characteristics.  To  employ  USVA  for  in  vivo  applications,  one  must  take  into  account  limitations 
such  as  safe  (ultrasound)  power  limit,  tissue  attenuation,  body  noise,  and  phase  aberration.  The 
required  ultrasound  power  to  generate  a  detectable  acoustic  emission  depends  on  the  object, 
acoustic  noise,  and  the  receiver  sensitivity.  Experimental  results  shown  in  Fig.  3  demonstrate 
that  ultrasound  intensities  as  low  as  30  mW/cm2  are  sufficient  to  detect  the  acoustic  emission 
from  a  1  mm  diameter  object  with  our  hydrophone  in  the  water  tank.  This  power  value  is  much 
smaller  than  the  FDA  limit  for  safe  diagnostic  ultrasound  applications.  Tissue  attenuation  reduces 
the  ultrasound  intensity  at  the  target,  and  hence,  the  acoustic  emission  (by  the  factor  A(zq)). 
Attenuation  limits  the  usable  ultrasound  frequency,  and  hence,  lowers  the  resolution.  It  also  limits 
the  signal-to-noise-ratio  (SNR)  because  of  the  loss  in  the  acoustic  emission  energy  as  a  result 
of  ultrasound  attenuation  by  tissue.  Direct  effect  of  tissue  attenuation  on  the  acoustic  emission 


signal  is  probably  negligible  because  attenuation  of  the  compressional  waves  at  frequencies  in 
the  order  of  a  few  kHz  in  soft  tissues  is  low.  Sources  of  biological  noise  of  the  human  body 
include  cardiovascular  and  respiratory  systems,  and  muscle  movements.  Body  noise  is  usually 
concentrated  below  1  KHz  and  can  be  filtered  out  if  A /  is  above  this  value.  The  SNR  can 
be  improved  by  increasing  the  time  duration  of  the  signal  recorded  at  each  point  (to  increase 
the  signal  energy)  and  using  very  narrow  band  filters  (to  reject  the  noise),  while  keeping  the 
ultrasound  power  within  the  safe  level.  Phase  aberration  in  tissue  can  reduce  the  sensitivity  of 
the  system  by  decreasing  the  effective  ultrasound  energy  density  at  the  beam  interaction  region. 
One  may  employ  known  phase  aberration  correcting  methods  to  reduce  such  an  effect.  The 
practical  value  of  these  methods  for  USVA  remains  to  be  studied. 

Material  evaluation — Another  field  in  which  USVA  can  be  potentially  useful  is  material 
characterization,  including  mechanical  parameter  evaluation,  imaging,  and  nondestructive  testing 
of  materials.  The  method  described  here  can  be  used  for  detection  and  imaging  flaws  in  materials. 
Also,  one  may  use  USVA  to  evaluate  the  mechanical  frequency  response  of  an  object  at  low 
frequencies.  In  such  case,  we  are  interested  in  determining  Qauj(x,v)  versus  frequency.  We 
assume  that  the  object  is  uniform  within  the  projected  area  S.  Then,  the  total  radiation  force  on 
this  object,  F^i,  can  be  calculated  by  integrating  F^u(x,y)  over  S.  Referring  to  Eqs.  21,  one 
can  show  that,  for  Aco  «  u0,  F&u  is  virtually  independent  of  Au/.  If  Ha W(Z)  is  known  and 
nonzero,  then  the  function  Qau(x,  y)  can  be  estimated  using  Eq.  16  as  Qau(x,  y)  =  ■J>*uJlx'v}n. 
In  practice,  Pau{x,v)  is  obtained  by  sweeping  Au>  in  the  range  of  interest  and  recording  the 
resulting  acoustic  emission  [12]. 
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9  LEGENDS 


Figure  1  Principle  of  ultrasound-stimulated  vibro-acoustography. 

Figure  2  Ultrasound-stimulated  vibro-acoustography  system.  The  confocal  ultrasound  annular 
array  transducer  with  two  elements  is  shown  on  the  left. 

Figure  3  Acoustic  emission  field  intensity  versus  the  combined  ultrasound  intensity. 

Figure  4  USVA  images  of  a  380  //.m  glass  bead:  (A)  in-phase,  (B)  quadrature,  (C)  phase,  and 
(D)  magnitude.  The  phase  in  (C)  ranges  from  -7r  radians  (black  regions)  to  +7r  radians  (white 
regions),  and  it  was  normalized  to  be  zero  at  the  center  of  the  glass  bead.  (Modified  with 
permission  from  [12].  Copyright  1998  American  Association  for  the  Advancement  of  Sciences.) 

Figure  5  The  theoretical  PSF  profile  of  the  USVA  system  according  to  Eq.  24  and  the  glass  bead 
in-phase  image  profile  (Fig.  4a)  obtained  from  the  experiment. 
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Acoustic  Field  Calculation  in  Ultrasound-Stimulated  Vibro-Acous- 
tography 
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Background:  In  ultrasound-stimulated  vibro-acoustography,  two  ultrasound 
beams  are  used  to  produce  oscillating  radiation  force  on  the  object  under  test. 
The  resulting  acoustic  field  is  used  to  image  the  object.  Here  we  discuss  the 
analytical  relationships  between  the  ultrasound  pressure  field  and  the  resulting 
acoustic  field  (acoustic  emission).  These  relationships,  which  involve  object 
mechanical  parameters,  can  be  used  to  deduce  object  characteristics  from  the 
received  acoustic  field. 

Method:  The  radiation  force  of  two  intersecting  ultrasound  beams  on  an 
object  is  calculated.  Assuming  linearity,  object  displacement  due  to  this 
force  and  the  resulting  acoustic  emission  pressure  field  are  calculated. 

Results:  It  is  shown  that  the  acoustic  pressure  amplitude  is  proportional  to 
input  ultrasound  power.  Also,  acoustic  amplitude  is  a  function  of  object 
mechanical  parameters  such  as  the  stiffness,  mass,  and  the  area. 
Conclusion:  The  ultrasound-stimulated  acoustic  emission  field  of  an  object 
is  related  to  object  mechanical  parameters  via  a  closed  form  relationship. 
This  allows  estimation  of  object  parameters  from  remote  measurements  of 
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Potential  Applications  of  Ultrasound- 
Stimulated  Vibro-Acoustography  in  Medical 
Diagnosis  and  material  Characterization 

Mostafa  Fatemi  and  James  F.  Greenleaf 
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Elastic  features  of  tissue  are  known  to  be  highly  correlated  with  tissue  state 
and  with  some  disease  processes.  Elastic  constants  are  linked  to  numerous  material 
properties,  and  hence,  are  valuable  in  nondestructive  evaluation  of  materials. 
Recently,  a  new  technique  for  imaging  <pbjects  and  evaluating  their  mechanical 
parameters  was  introduced  by  the  authors  [Science,  (280):82-85,  1998].  This 
technique,  ultrasound-stimulated  vibro-acoustography  (USVA),  uses  ultrasound 
to  produce  an  oscillating,  highly  localized  radiation  force  to  vibrate  the  object 
at  low  frequencies.  The  resulting  acoustic  emission  is  a  function  of  the  local 
ultrasonic  properties,  which  are  measured  by  the  conventional  ultrasound,  and 
are  also  a  function  of  mechanical  parameters  (such  as  the  elastic  constants).  This 
method,  therefore,  emulates  the  simplicity  and  high  resolution  of  the  conventional 
ultrasound,  while  providing  additional  information  on  low  frequency  mechanical 
features  of  the  object.  In  the  medical  field,  this  method  can  provide  high  resolution 
images  of  tissue  elasticity.  Estimation  of  elastic  parameters,  such  as  stiffness,  are 
valuable  in  diagnosis,  for  instance,  to  detect  the  presence  of  calcification  within 
arteries  or  breast.  Remote  and  nondestructive  evaluation  of  viscosity,  resonant 
frequency,  and  dynamic  Young’s  modulus  are  potential  applications  of  USVA  in 
material  science.  Experimental  results  demonstrating  various  applications  will  be 
presented. 
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Beam  Forming  for  Ultrasound-Stimulated  Yibro-Asoustography 
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Ultrasound-Stimulated  Vibro-Asoustography  (USVA)  is  a  recently  developed  method  that 
utilizes  ultrasound  energy  to  evaluate  the  low-frequency  mechanical  parameters  of  an  object 
[Science  (280):82-85,  1998],  The  novelty  in  this  method  results  from  combining  the  benefits 
of  ultrasound  with  those  of  low-frequency  elastic  waves.  This  method  takes  advantage  of  the 
radiation  force  phenomenon  to  directly  convert  ultrasound  energy  to  low  frequency  elastic 
waves  in  the  object.  Hence,  the  object  response,  in  general,  is  a  function  of  both  the 
elastic  parameters  and  ultrasonic  characteristics  of  the  object.  The  low  frequency  features 
of  the  object  is  of  particular  interest  in  many  applications,  including  tissue  imaging  and 
characterization,  and  material  science.  Elastic  parameters  such  as  object  stiffness,  which 
are  not  available  from  the  conventional  ultrasonic  methods,  are  valuable  in  medical  field,  for 
example  to  delineate  hard  deposits  of  calcifications  in  arteries  from  the  surrounding  soft  tissue. 

The  principles  used  for  ultrasound  beam  forming  in  USVA  is  fundamentally  different 
from  those  for  the  conventional  ultrasound  imaging  methods.  The  reason  is  that  in  USVA  the 
ultrasound  field  is  designed  to  interact  nonlinearly  with  the  object  whereas  in  the  conventional 
ultrasound  they  interact  linearly. 

In  this  paper,  we  describe  USVA  beam  forming  design  methods,  and  present  three  beam 
forming  configurations,  including  single,  dual  cross-beam  (x-focal),  and  confocal  beams.  In 
the  latter  two  configurations,  the  two  ultrasound  sources  are  driven  at  different  frequencies. 
The  acoustic  field  pressure  and  the  resulting  radiation  stress  on  an  elementary  object  are 
calculated.  Analytical  and  simulation  results  indicate  that  the  pressure  and  stress  fields 
produced  by  the  x-focal  and  confocal  beams  are  circularly  asymmetric  and  symmetric, 
respectively.  The  stress  field  in  the  confocal  case  is  distributed  according  to  the  product 
of  two  Bessel  functions  of  the  first  kind.  Simulation  results  indicate  that,  given  transducers 
with  equal  aperture  size,  the  confocal  beam  produces  the  narrowest  central  lobe  while  the 
x-focal  beam  has  the  lowest  sidelobe  amplitude.  Experimental  results  using  focused  3.0  MHz 
transducers  in  a  water  tank  are  in  agreement  with  the  theoretical  and  simulation  findings. 


